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The  work  described  in  this  report  was  conducted  by  the 
McDonnell  Douglas  Astronautics  Company  under  Independent 
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lished  as  an  Air  Force  technical  document  because  the  material 
represents  the  analytical  framework  that  is  a  central  part  of  the 
"Rocket  Propulsion  System  Hardening  Design  Handbook,  Vol¬ 
ume  1,  Nuclear"  developed  by  McDonnell  Douglas  Astronautics 
Company  for  the  Air  Force  Rocket  Propulsion  Laboratory,  under 
Contract  No.  F046  1 1  - 73 -C -0042. 
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Beale . 
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1  g/cmJ  (units  of  weight  density)  =  1  dyne -sec  /cm  (units  of  mass  density) 
1  gram  =  1  dyne-sec^/cm  (=  980.  7  dynes  at  acceleration  of  1  gravity) 

1  dyne  =  1  gram  -  cm/sec^ 

1  erg  =  1  dyne  -  cm  =  10"^  joules 

1  calorie  =  4.  186  x  10^  ergs  =  4.  186  x  10"^  megabar-cm^ 

1  tap  =  1  dyne-sec/cm^  =  1  bar-psec 

1  bar  =  10&  dynes/cm^  =  14.  5  psi  =  0.  987  atmosphere 

1  psi  =  69,  000  dynes/ cm2  =  0.69x  10_?  megabar  =  70.  3  1  g/cm^ 

1  kilobar  =  10^  bars  =  10^  dynes/cm^ 

1  megabar  =  10-5  kilobars  =  10  dynes/cm 
1  cal/ g°C  =  1  Btu/lb°F 

1  ft/sec  =  30.  48  cm/sec  =  30.48  x  10-&  cm/psec 
1  kev  =  10^  electron  volts  =  1.  603  x  10-9  ergs 
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Section  1 

INTRODUCTION  AND  SUMMARY 

This  report  summarizes  analytical  relationships  derived  for  estimating 
the  effects  of  x-rays  on  materials.  Since  all  major  x-ray  effects  are  treated 
in  a  general  way,  they  are  generally  applicable  to  anything  exposed  to 
x-rays.  Although  the  relationships  were  developed  primarily  for  the  designer 
and  systems  analyst,  their  quality  is  comparable  to  computer  codes  so  they 
are  also  useful  to  the  technologist. 

The  effects  covered  include  energy  deposition  in  single  and  multiple 
materials,  front-face  mass  removal,  structural  heating,  blowoff  impulse, 
stress  generation  by  impulse  loads  and  thermal  loads,  stress  wave  attenua¬ 
tion  in  undistended  and  distended  materials,  and  structural  response.  Opti¬ 
mum  energy  spectra  which  maximize  an  effect  for  a  given  fluence  are 
emphasized;  however,  specific  spectra  are  also  treated. 

Except  for  stress  wave  attenuation  in  distended  materials,  the  relationships 
and  their  demonstrated  general  validity  for  effects  analysis  are  new — 
they  were  developed  during  the  past  year.  The  motivation  for  developing  an 
analytical  approach  was  to  circumvent  the  difficulties  associated  with  : 

(1)  the  need  to  provide  practical  understanding  of  individual  and  combined 
x-ray  effects  to  individuals  not  conversant  with  x-ray  effects,  (2)  the  numerous 
factors  involved  in  analyzing  x-ray  effects,  (3)  the  vast  number  of  materials 
and  material  combinations  associated  with  some  systems,  (4)  the  determination 
of  spectrally  optimized  effects,  (5)  the  need  for  flexibility  in  accommodating 
refinements  in  damage  or  failure  criteria,  (6)  the  tedium  and  crudeness  of 
existing  analytical  techniques,  (7)  the  designer's  lack  of  convenient  access 
to  computerized  analytical  techniques,  and  (8)  the  inherent  problems  of  over¬ 
dependence  on  computerized  numerical  analysis. 
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The  resulting  generalized  analytical  relationships  are  simple  closed-form 
equations.  With  only  three  exceptions,  they  are  all  derived  from  basic 
theoretical  concepts  rather  than  empirical  correlations;  otherwise,  the 
task  would  have  been  impractical.  They  have  been  compared  with 
computer  analyses  performed  with  energy-deposition  codes  and  hydro- 
dynamic  codes,  and  the  comparisons  have  shown  the  equations  to  produce, 
over  the  complete  practical  ranges  of  parameters,  the  results  of  the  computer 
codes  to  within  generally  10  to  30  percent.  It  is  not  clear  how  much  of  this 
difference  is  due  to  approximations  in  the  equations  and  how  much  is  due  to 
errors  in  the  codes.  In  some  cases,  the  error  appears  to  be  in  the  codes, 
for  example,  in  consistent  discontinuities  in  computer  results  which  have 
no  physical  basis,  and  the  error  becomes  apparent  only  when  a  parametric 
correlation  is  attempted.  When  the  differences  between  calculations  that 
involve  a  number  of  parameters  are  small,  it  takes  considerable  time  and 
effort  to  determine  the  exact  sources  of  the  differences.  There  have  been  a 
few  isolated  cases  when  the  difference  between  analytical  results  and  com¬ 
puter  results  has  been  as  much  as  a  factor  of  2  (e.  g.  ,  the  empirical  energy 
deposition  relationships);  however,  for  the  same  or  similar  cases,  two  com¬ 
puter  programs  gave  results  that  differed  by  greater  than  a  factor  of  2,  pre¬ 
sumably  because  they  use  different  sets  of  x-ray  mass  absorption 
coefficients. 

In  adopting  a  perspective  for  the  accuracy  of  any  analytical  tool  for 
vulnerability  analysis,  it  is  helpful  to  realize  that,  from  a  practical  stand¬ 
point,  useful  precision  is  limited  by  the  accuracy  of  materials  data,  failure 
criteria,  and  the  reality  of  assumptions  regarding  radiation  environments. 

The  uncertainties  of  all  of  these  far  exceed  the  uncertainty  in  the  present 
analytical  modeling  approximations. 

A  large  number  of  comparisons  have  been  made.  However,  because  of 
the  numerous  relationships  and  ranges  of  parameters  involved,  it  has  not 
been  practical  to  make  a  comprehensive  study  of  the  accuracy  of  the  analyti¬ 
cal  relationships.  It  has  been  necessary  to  present  the  results  of  the  com¬ 
parisons  that  have  been  made  in  somewhat  general  terms;  documentation  has 
not  been  possible  due  to  the  limited  scope  of  effort.  Nevertheless, 
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sufficient  comparisons  have  been  made  over  the  full  ranges  of  parameters 
of  the  individual  relationships  that  their  validity  and  nominal  accuracy  have 
been  established. 

The  objective  has  been  to  duplicate  computer  analyses  as  closely  as  possible. 
Aside  from  the  benefits  of  an  uncomplicated  analysis,  such  an  approach  gives 
numerical  results  which  have  the  same  state-of-the-art  limitations  relating 
to  material  properties,  failure  criteria,  and  theoretical  constructs  as  the 
computer  programs  it  seeks  to  emulate.  These  same  limitations  on  the 
application  of  theoretical  tools  make  them  indispensable  to  a  rational  design 
approach  because  of  the  even  more  severe  restrictions  they  impose  on  the 
alternate  approach  of  empiricism^. 

In  addition  to  presenting  the  derivation  of  equations  for  each  phenomenon, 
the  last  section  of  the  report  combines  these  equations  to  obtain  explicit 
expressions  for  fluence  thresholds  for  various  failure  mechanisms  and 
situations  relevant  to  propulsion  systems.  Others  can  be  similarly  derived. 

An  attempt  has  been  made  to  maintain  consistent  units  and  symbols 
throughout.  They  are  defined  in  each  context.  Units  of  calories,  grams, 
centimeters,  microseconds,  kilobars,  and  kilotaps  have  beer  used  consist¬ 
ently.  A  set  of  unit  conversion  relations  has  been  included  for  reference. 


J.  Watcher,  Rocket  Propulsion  System  Hardening  Design  Handbook. 
McDonnell  Douglas  Paper  WD  23  53.  Presented  at  the  Fourth  Symposium 
on  Nuclear  Survivability  of  Propulsion  and  Ordnance  Systems,  Defense 
Nuclear  Agency  Report,  1974. 
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Section  2 


X-RAY  ENERGY  DEPOSITION 


The  most  fundamental  quantity  required  for  x-ray  effects  analysis  is  the 
x-ray  energy  deposition  in  a  given  material  or  laminate;  the  analysis  of  all 
effects  relies  upon  this  information.  An  x-ray  spectrum  is,  for  convenience, 
often  characterized  by  a  single  blackbody  temperature  expressed  in  kev 
(10  electron  volts).  An  x-ray  environment  can  also  be  characterized  by  a 
range  of  blackbody  temperatures.  The  objective  of  this  section  is  to  derive 
useful  analytic  relationships  for  energy  deposition  for  blackbody  energy 
spectra. 

2.  1  MAXIMUM  DOSE  FOR  A  RANGE  OF  BLACKBODY  SPECTRA 
2.  1.  1  Single  Material 

For  a  range  of  blackbody  temperatures  and  for  a  given  point  within  a 
material,  a  unique  quantity  is  the  maximum  dose  (i.  e.  ,  energy  depostion)  for 
a  given  fluence.  The  exact  expression  for  the  dose  in  a  single  material  is 

D(x,  Z,  T)  =  4>  /  B(E,  T)  exp  [  -  p(Z,  E)p  x]  p.  (Z,  E)  dE  (1) 

Jo 


where  D  is  the  dose  in  cal/g,  x  is  the  thickness  variable  in  cm,  Z  is  the 

atomic  number  of  the  material,  T  is  the  blackbody  temperature  in  kev, 

?  3 

4>  is  the  incident  x-ray  fluence  in  cal/cm  ,  p  is  the  material  density  in  g/cm  , 

2 

p(Z,  E)  is  the  x-ray  mass  absorption  coefficient  in  cm  /g,  E  is  the  x-ray 
(photon)  energy  in  kev,  and  B(E,  T)  is  the  blackbody  function  given  by 


B(E,  T)dE 


(2) 


where  y  -  E/T.  The  integral  of  Equation  (2)  from  zero  to  infinity  is  one. 


Disregarding  absorption  edges,  it  is  well  known  that  a  good  approximation 
for  the  mass  absorption  coefficient  is 

p(Z,  E)  -  ZnE"m  (3) 

where  jjl ^  is  the  reference  value  at  1  kev.  For  the  range  of  energies  of 
practical  interest,  n  «  m  «  3. 

Defining  the  variable  =  kj  Zn  E  m  px,  Equation  (1)  becomes 


The  objective  is  to  maximize  D  by  maximizing  the  integral.  Since  e  ^  dll' 
varies  far  more  strongly  with  y  than  does  the  blackbody  function  in  paren¬ 
theses,  a  good  assumption  is  to  remove  the  blackbody  function  from  the 
integral  and  maximize  it,  with  the  remainder  of  the  integral  becoming  unity. 
Thus, 


Max 


,4  -4 
4  e 


(5) 


Substituting  this  and  m  =  3  into  Equation  (4)  gives  for  the  maximum  dose  for 
any  fluence,  in  any  material, 


D 


m 


4  -4 

1  5  (4)  e 


3tt 


Jl  = 

px 


0.  24 

px 


(6) 


Because  this  is  an  optimum,  the  variation  of  dose  with  all  parameters  is 
weakest  near  D  .  Also,  the  floating  parameter  which  is  optimized  (i.  e.  ,  y) 
has  the  effect  of  obscuring  the  absorption  edges  in  higher- atomic-number 
materials  so  that  they  are  manifested  only  in  terms  of  an  average  m,  a 
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parameter  that  is  a  weak  function  of  Z  and  a  weaker  function  of  E.  This  is 

verified  by  parametric  computer  calculations  with  a  typical  one- dimensional 

point-kernel  (exponential  attenuation)  x-ray  transport  code  using  the  mass 

2 

absorption  coefficients  of  Marotta. 

Table  1  lists  coefficients  of  Equation  (6)  that  were  determined  by  parametric 
computer  calculations.  The  coefficients  naturally  fall  into  three  categories: 
Group  1,  Z  s  6;  Group  2,  12  S  Z  S  30;  and  Group  3,  Z  2  40.  However,  the 
single  coefficient  of  0.  24  is  within  ±20  percent  for  all  materials,  and  is 
sufficiently  accurate  for  virtually  all  situations  which  do  not  fall  outside  a 
specified  range  of  T.  (What  to  do  when  they  do  fall  outside  is  discussed  in 
Section  2.  2).  Equation  (6)  is  therefore  applicable  to  any  homogeneous 
material,  and  is  independent  of  all  physical  properties  except  density. 

Greater  accuracy  is  obtained  by  using  the  coefficients  of  Table  1,  but  it  is 
helpful  to  realize  that  since  blackbody  spectra  or  any  other  idealized  spectra 
represent  convenience  more  than  reality,  consistency  is  more  important  than 
absolute  accuracy  once  a  practically  acceptable  level  of  accuracy  is  attained. 

The  variation  of  the  mass  absorption  coefficient  is  shown  in  Figure  1  for 
carbon,  aluminum,  iron,  and  tungsten.  The  sharp  discontinuities  caused  by 
the  various  or bital- electron  shells  (designated  K,  L,  M,  etc.)  are,  in 
effect,  completely  smoothed  out  by  optimizing  D  in  terms  of  a  dispersive 
spectrum.  For  energies  higher  than  the  K-shell  potential,  m  =  3.  0  for  all 
elements,  and  for  all  lower  energies,  m  =  2.  5.  Between  carbon  and  alum¬ 
inum,  n  =  3.4.  Between  aluminum  and  iron,  n  =  3  3.  And  between  iron  and 
tungsten,  n  =  3.  1  .  These  values  of  n  relate  only  to  energies  higher  than 
the  K-shell  potentials. 

Figure  1  indicates  several  things  relevant  to  x-ray  dose.  First,  carbon  is 
the  highest  -Z  material  of  interest  which  has  no  significant  absorption  edge 
(i.  e.  ,  for  E  1  kev);  thus,  other  materials  will  behave  somewhat  differently 
than  carbon  and  beryllium.  Second,  for  a  dispersive  spectrum,  the  effect 
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C.  R.  Marotta,  Updated  Master  Library  Tape  for  PHOTRAN,  Air  Force 
Weapons  Laboratory,  AFW  L-  TR-  67-  1 1 ,  1968. 


Table  1 


COMPUTER- CALCULATED  COEFFICIENTS  FOR 
MAXIMUM- DOSE  EQUATION 


i 

Element 

Atomic  Number 

Coefficient 

Beryllium 

.  . 

4 

0.  20 

Carbon 

6 

0.  20 

Aluminum 

13 

0.  22 

Titanium 

22 

0.  22 

Iron 

26 

0.  22 

Nickel 

28 

0.  22 

Copper 

29 

0.  22 

Niobium 

41 

0.  29 

Molybdenum 

42 

0.  29 

Tungsten 

74 

0.  29 

CR7S 


tr4  io°  101  io2 

E  (KCV) 

Figure  1.  X-Ray  Mass  Absorption  Coefficients 
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of  absorption  edges  is  minimal  except  for  the  average  absorption  coefficient 
characteristics  in  terms  of  effective  values  of  m,  and  n.  Third, 
dose  correlations  based  on  the  behavior  of  the  absorption  coefficient  at 
energies  higher  than  the  K- shell  potential  are  not  valid  for  Z  ~  12  because 
of  the  significant  influence  of  lower  energies;  the  actual  dependence  on  Z 
and  E  is  weaker  than  for  m  «  n  «  3.  Fourth,  the  effective  absorption 
coefficients  have  similar  behavior  only  within  the  three  categories  of  atomic 
number  mentioned  above,  and  in  the  highest  one  the  effective  absorption 
coefficient  is  approximately  the  same  for  all  materials,  with  the  main 
difference  being  the  individual  shell  potentials,  while  their  collective  effect 
is  nearly  identical.  Finally,  an  important  property  of  x-ray  absorption 
coefficients  is  that  they  vary  smoothly  with  atomic  number,  including  shell 
potentials,  and  therefore  integral  quantities  like  doses  for  dispersive  spectra 
also  vary  smoothly  with  atomic  number. 

Figures  2  through  7  show  dose  profiles  for  unit  fluence  and  for  1  to  1  5  kev 
blackbody  spectra  in  beryllium,  carbon,  aluminum,  iron,  niobium,  and 
tungsten.  The  maximum  dose  per  fluence  at  a  given  depth  is  represented  by 
the  envelope  formed  by  these  curves  (which  is  more  evident  for  the  higher  Z 
materials).  Also  shown  in  the  figures  are  plots  of  the  cumulative  fraction  of 
fluence  absorbed  to  a  given  depth.  These  data  were  generated  with  the 
previously  mentioned  computer  program.  The  curves  can  be  interpolated  in 
terms  of  the  effective  atomic  number  (see  Equation  (10))  and  px  (g/cm  )  to 
handle  any  homogeneous  material. 

It  is  important  to  realize  that  these  computer-generated  data  and  any  similar 

data  do  not  represent  absolute  accuracy.  For  example,  the  energy  deposition 

routine  in  the  PUFF  code,  which  uses  the  mass  absorption  coefficients  of 

3 

the  National  Bureau  of  Standards  ,  gives  results  that  in  some  cases  differ  by 
more  than  a  factor  of  2  from  the  data  in  the  figures.  This  is  presumed  to  be 
due  to  a  difference  in  the  absorption  coefficients. 


^G,  R.  White,  National  Bureau  of  Standards,  Report  NBS  1003,  1952. 
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Figure  2.  Energy  Deposition  in  Beryllium  for  Blackbody  Spectra 
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Figure  3.  Energy  Deposition  in  Carbon  for  Blackbody  Spectra 
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Figure  4.  Energy  Deposition  in  Aluminum  for  Blackbody  Spectra 
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Figure  5.  Energy  Deposition  in  Iron  for  Blackbody  Spectra 


CR79 


DEPTH (CM) 


Figure  6.  Energy  Deposition  in  Niobium  for  Blackbody  Spectra 
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f.  CUMULATIVE  DOSE  [(CAL/CM2)/(CAL/CM2)] 
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2.  1.  2  Laminated  Materials 


For  an  absorber  material  of  atomic  number  Z  and  any  intervening  material 
(generally  called  shielding)  of  atomic  number  Z  ,  which  have  the  same 
effective  values  of  m  and  n.  Equation  (1  )  Leads  directly  to  Equation  (6) 
except  for  modification  by  the  factor  (Zn/Zn)  to  account  for  the  scaled 
absorption  coefficient  as  defined  by  Equation  (3). 

That  is, 


D  =  0.  24  -$-(zn/Zn) 
m  Pk  \  a  s) 


(7) 


Parametric  computer  calculations  for  Z  <30  and  Z  <  30,  and  with  n  =  3.  3, 

3  S 

compare  with  Equation  (7)  to  within  better  than  20  percent. 


Because  for  Z  >  40  the  effective  absorption  coefficient  is  approximately  the 
same  for  all  materials  for  the  energy  range  of  interest  (see  Figure  1), 
parametric  computer  calculations  show  that  is  independent  of  material 
for  Z  >  40.  That  is,  Equation  (6)  is  also  applicable  for  any  combination  of 
materials  for  Z  >  40  to  within  better  than  20  percent  for  the  materials 
checked  (ranging  from  niobium  to  tungsten). 


Many  parametric  computer  calculations  were  made  of  numerous  group  3 
materials  over  numerous  group  1  and  group  2  materials.  Because  the 
effective  m's  and  n's  are  different  for  these  two  categories  of  materials, 
energy- dependence  and  material-dependence  are  no  longer  separable,  and 
consequently  the  functional  dependence  will  be  different  from  the  above 
relationships,  requiring  an  empirical  correlation.  It  was  found  that  for 
Za  <  30  and  Zg  2  40,  a  simple  empirical  expression  for  is 

D  =  0.40  - *  (zn/Zn\ 

m  (px)3/2  l  a  s  ' 

where  n  =  2.  8  and  the  correlation  is  within  *40  percent 
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Table  2  summarizes  the  parameters  for  the  general  expression 

Dm  =  a*  (pxfb  (z“/z“)  (9) 

where  a,  b,  and  n  are  fixed  for  the  material  combinations  defined  in  the 

table-  Parameters  have  not  been  determined  for  Z  >  40  and  Z  <  30. 

a  s 

2.  1.  3  Effective  Atomic  Number 

The  above  discussion  has  been  in  terms  of  pure  elemental  materials.  How¬ 
ever,  the  relationships  have  been  shown  by  computer  calculations  to  hold 
for  alloys  and  composite  materials  as  well.  In  fact,  the  above  expressions 

hold  for  any  material  or  combination  of  materials  if  Z,  Z  ,  and  Z  are  each 

3.  S 

represented  by  an  effective  atomic  number  defined  as 


(10) 


where  w^  is  the  mass  fraction  of  each  element,  and  the  sum  is  over  all 

elements  in  either  the  shield  or  the  absorber.  Whether  the  shield  is  a 

homogeneous  material  or  a  laminate  is  irrelevant.  The  shield  is  defined  as 

2 

the  total  mass  of  material  (in  g/cm  )  between  the  incident  radiation  and  the 
point  at  which  the  dose  is  evaluated  (including  air).  The  value  of  n  =  3  in 
Equation  (10)  adequately  represents  all  material  combinations  for  the  purpose 
of  determining  Z  ;  very  nearly  the  same  value  of  Z  is  obtained  for 
2.  8  <  n  —  3.  3.  Because  Z  is  a  strong  function,  Z^  is  usually  dominated  by 
one  or  two  elements  in  a  material. 


If  a  laminated  "shield"  is  composed  of  M  different  materials  whose 

individual  effective  atomic  numbers  are  known,  a  modified  version  of 

Equation  (10)  is  convenient  for  determining  the  effective  atomic  number  of 

the  shield,  3,  . 

**» 


\!zVi 


2  1 


r  3 


1/3 


+  PMXMZ  eM 


‘  2 


'  MXM 


(11) 
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Table  2 

PARAMETERS  FOR  MAXIMUM  -  DOSE  EQUATION* 


Z  Range 

For  Shield 

Z  Range 

For  Absorber 

Coefficient 

a 

Z  Exponent 
n 

Thicknes  s 
Exponent 
b 

Z  <  30 

Z  <  30 

0.  24 

3.  3 

1 

Z  >  40 

Z  >  40 

0.  24 

0 

1 

Z  >  40 

Z  <  30 

0.  40 

2.  8 

3/2 

Z  <  30 

Z  >  40 

? 

? 

? 

Any  Z 

Z  =Z 
a  s 

0.  24 

- 

1 

♦General  Equation:  =  a4>  (p 

x)-b  (z“/z») 

Figure  8  is  convenient  for  relating  the  Z  factor 

number  (Z  or  Z  ). 

e 


or  Z 


:) 


to  the  atomic 


101  102  103  104  10®  106 

z" 


Figure  8.  Plot  of  Zn  versus  Z 
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2.  2  DOSE  FOR  SPECIFIC  SPECTRA 

For  any  specific  spectrum,  an  effective  general  method  of  dose  analysis  is 

based  on  the  above  analytical  approach.  However,  it  is  not  completely 

analytical,  but  requires  a  number  of  dose-versus-depth  (in  g/cm  )  curves 

for  the  range  of  effective  atomic  numbers  of  interest.  If  specific  blackbody 

spectra  are  of  interest,  the  parametric  curves  are  like  those  in  Figures  2 

through  7,  except  it  is  more  convenient  if  the  material  depth  is  expressed  in 
2 

g/cm  instead  of  cm. 

A  conveniently  small  number  of  curves  can  represent  the  full  range  of 
and  T  for  blackbody  spectra;  that  is,  with  up  to  20  to  30  figures  like 
Figures  2  to  7,  a  single  homogeneous  material  with  a  particular  Zg  can  be 
adequately  represented  by  the  parametric  energy-deposition  curves  of  the 
material  with  the  closest  Z^,  without  interpolation  in  Z^  being  required.  It 
is  important  to  represent  one  material  by  another  in  terms  of  g/cm^  and  not 
in  cm  to  account  for  differences  in  density. 

For  laminated  materials,  the  same  procedure  is  applied  to  determine  the 
"shield  dose,  "  Ds,  in  the  homogenized  representation  of  the  shield  with 
effective  atomic  number  Z  .  The  dose,  D  ,  in  the  absorbing  material  behind 

6S  cl 

the  shield  is  then  obtained  by  scaling  D  by  the  ratio  of  Z  factors: 

s 


Thus,  with  a  set  of  parametric  energy- deposition  curves  for  enough  single 
materials,  any  material  configuration  can  be  easily  analyzed.  The  appropri¬ 
ate  values  of  n  for  Equation  (12)  are  obtained  from  Table  2. 

2.  3  EXPONENTIAL  REPRESENTATION  OF  DOSE  PROFILES 
2.  3.  1  Basic  Exponential  Function 

Most  actual  dose  profiles,  including  those  of  blackbody  spectra,  exhibit 
exponential- like  behavior  down  to  a  significant  degree  of  dose  attenuation. 

To  the  extent  that  an  exponential  function  adequately  represents  the  portion 
of  the  dose  profile  that  is  dominant  for  a  given  effect,  closed- form  analytical 
solutions  can  be  obtained  for  estimating  most  x-  ray  effects.  This  approach 
was  investigated  and  found  to  be  successful.  For  optimized  quantities 
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involving  the  integral  of  the  dose  (e.  g.  ,  blowoff  impulse),  completely 
satisfactory  results  can  be  obtained  even  when  an  exponential  function  is  a 
poor  approximation  of  the  dose  profile. 

The  general  representation  of  dose  profiles  as  an  exponential  function  is 

jc 

D  =  D  e  *  (13) 

o 


where  D  is  the  surface  dose  and  £  is  the  characteristic  dose  attenuation 
o 

depth  in  cm,  or  the  inverse  of  the  effective  (average)  linear  attenuation 
coefficient  p  (i.  e.  ,  l  -  p  /P-). 


To  conserve  energy,  the  integral  of  Equation  (13)  over  all  x  must  equal  the 
incident  fluence.  This  holds  only  for  an  exact  exponential  profile  and  results 
in  a  unique  £  which  is  related  to  by 


£ 

o 


4> 


(14) 


However,  in  some  cases  (e.  g.  ,  thermomechanical  stress)  the  quantity  to  be 

approximated  is  not  the  total  energy  absorbed,  but  the  initial  part  of  the 

dose  profile.  In  this  case,  £  is  defined  as  the  thickness  required  to  reduce 

the  dose  to  e  *D  . 

o 


Figure  9  shows  how  the  dose  profiles  for  1,  5,  and  1  5  kev  typically  vary 
when  divided  by  Dq  and  plotted  in  terms  of  x/f^.  Also  shown  are  the  two 
kinds  of  exponential  repr esentations  for  £  and  £  ,  with  the  areas  under  the 
curves  indicating  energy  is  not  conserved  for  £*£  .  It  is  seen  that 
exp  (-x/i)  approximates  the  actual  profile  well  down  to  D/D^  =  0.  Z,  and 
underestimates  the  dose  by  about  a  factor  of  2  at  D/Dq  =0  1,  Comparisons 
with  computer  calculations  indicate  that  this  agreement  to  this  degree  of 
attenuation  is  adequate  for  predicting  x-ray  effects  for  most  situations  if 
conservation  of  energy  is  introduced  when  it  is  important. 
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t  o 

Figure  9.  Comparison  of  Blackbody  Dose  Profiles  with  Exponential  Approximation 

2.  3.  2  Exponential  Parameters  for  Blackbody  Spectra 

To  apply  Equation  (13)  and  relationships  which  are  derived  from  it  to 

specific  blackbody  spectra  or  any  other  spectra,  Dq  and  2  can  be  obtained 

from  parametric  data  like  those  shown  in  Figures  2  through  7.  However,  it 

is  useful  to  have  an  analytical  correlation  for  D  and  2.  Parametric  com- 

o 

puter  calculations  were  used  to  determine  such  expressions  empirically  for 
blackbody  spectra. 

The  surface  dose  is  defined  by  Equation  (1)  with  x  =  0.  This  gives 


D  = 
o 


4>  j  B(E,  T)p(Z,  E)dE  = 

J  r\ 


B(y)p1  dy/y 


m 


(15) 


For  Z  2  6,  |i^  is  a  constant,  m  3,  and  the  integral  is  a  constant.  However, 
for  Z  s  12,  this  is  no  longer  true  (see  Figure  1),  and  the  integral  is  a  function 
of  Z  and  T  which  must  be  determined  empirically.  The  following  relation¬ 
ships  were  found  to  correlate  all  computer  -  calculated  results. 


-  3  3 

3T  '  Z  <|>, 


Z  <  6 


(16a) 


D  =  60  Z  T  d>,  Z  >  12 
o  T 


(16b) 


Similarly,  the  correlations  for  &  are 


£  =  0.015 


Z  <  6 


(17a) 


3/2  /  ^  rl/4 


£  =  4X10' 


=*  4  X  10 


-4  T' 


Z  >  12 


(17b) 
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In  these  equations,  T  is  in  kev,  <|>  in  cal/cm  ,  p  in  g/cm  ,  Dq  in  cal/g,  and 
SL  in  cm. 


These  correlations  were  obtained  with  many  claculations  using  the  mass 

2 

attenuation  coefficients  of  Marotta  .  Equations  (16a)  and  (17a)  reproduce 
the  results  to  within  a  few  percent.  Equations  (1  6b)  and  (17b)  reproduce 
the  results  generally  to  within  30  percent.  However,  there  are  a  few  dis¬ 
continuities  in  the  calculated  data  in  addition  to  apparently  random  scatter, 
which  disagree  with  the  equations  by  as  much  as  a  factor  of  2.  This  is 
attributed  to  uncertainties  in  the  mass  absorption  coefficients,  or  the  numer¬ 
ics  of  the  calculation,  or  both,  since  D  and  i  must  both  be  smooth  functions 

o 

of  Z  and  E  for  dispersive  spectra.  (Several  comparisons  of  Equations  (16) 

and  (17)  with  results  obtained  with  a  version  of  the  PUFF  computer  program, 

which  uses  the  mass  absorption  coefficients  of  the  National  Bureau  of  Stand- 
3 

ards  ,  indicate  somewhat  better  agreement  than  with  the  Marotta  data,  but 
indicate  the  same  general  agreement.  However,  in  a  few  cases  the  compari¬ 
son  is  not  as  good.  It  is  not  known  which  computer  results  are  best.) 
Because  of  the  scatter  in  the  calculated  data,  the  approximation  of  Equa¬ 
tion  (17)  which  neglects  Z- dependence  has  the  same  relative  accuracy  as 
when  the  Z-dependence  is  retained. 


'Marotta,  op.  cit. 
1  White,  op.  cit. 
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2.  4  OPTIMUM  TEMPERATURE  FOR  MAXIMUM  DOSE 

Because  there  are  always  practical  limits  to  the  range  of  blackbody  tempera¬ 
tures  to  be  considered,  it  is  necessary  to  know  the  optimum  temperature 
corresponding  to  a  particular  Dm  in  order  to  determine  if  it  is  within  the 
range  of  interest.  If  it  falls  outside  the  range,  then  the  dose  and  associated 
quantities  are  estimated  on  the  basis  of  the  appropriate  limiting  temperature. 


To  obtain  the  optimum  temperature  T 


op 


is  derived  for  an  exponential 


dose  profile.  This  is  done  by  equating  to  zero  the  differential  of  Equation  (13) 

with  respect  to  £  ,  using  Equation  (14)  to  relate  i  and  D  ,  This  results  in 
o  o  o 

the  optimum  £  for  an  exact  exponential  profile: 


op 


(18) 


and 


D  =  - — £  =  0.  37  — 
m  p  x  p  x 


(19) 


It  is  first  noted  that  Equation  (18)  results  in  Equation  (19),  which  is  very 

similar  to  Equation  (6).  However,  it  cannot  be  used  directly  by  substituting 

Equation  (17)  for  £  and  solving  for  T  ,  nor  does  it  imply  that  D  /D  =  e, 
M  op  op  r  !  o  m 

both  of  which  are  far  from  correct.  Instead,  Equation  (18)  indicates  a 
general  integral  relationship  for  D^;  namely,  that  for  any  exponential- like 
profile,  the  cumulative  fractional  energy  deposition  f,  to  the  depth  at  which 
is  evaluated,  is  f  =  1  -e  ^  =  0.  63.  This  compares  with  computer  cal¬ 

culations  for  blackbody  spectra  to  within  10  percent.  If  f  (e  ^)  is  approxi¬ 
mately  constant  for  a  certain  class  of  dose  profiles,  then  so  is  the  total 
energy  balance  approximately  a  constant.  Consequently,  integrating 
Equation  (13)  results  in 


D  l 
o 


D  £ 
oe  o 


D  4> 
oe 

pD 

o 


(20) 
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where  D  is  the  surface  dose  corresponding  to  an  exponential  relationship, 

which  in  this  case  is  Equation  (19),  and  £  is  related  to  D  by  Equation  (14). 

Equation  (20)  is  the  general  expression  for  equating  the  total  energy  of  an. 

exponential  approximation  represented  by  Dq  and  £  to  that  for  an  exponentially 

derived  relationship  expressed  in  terms  of  D  .In  this  case,  D  is  given  bv 

oe  oe  ’ 


D  =  eD  =  -&- 
oe  m  p  x 


Substituting  this  into  Equation  (20)  results  in 


d2£  =  -i- 

o  2 

P  x 


(21) 


(22) 


Substituting  Equations  (16)  and  (17)  into  Equation  (22)  gives 


J/3 


Top  =  0.5  Z  (p x)  '  ,  Z  56 


Top  =  1. 2  Z4/5  (px)2/5,  Z  >  12 


(23a) 


(23b) 


Equation  (23)  is  valid  for  all  maximum- dose  relationships  for  which  b  =  1 

(see  Table  2).  The  atomic  number  of  the  absorber  does  not  influence  T  ; 

op 

the  Z  in  Equation  (23)  refers  to  the  effective  Z  of  the  shield  material, 
i.  e.  ,  Zg  as  defined  by  Equation  (11). 

For  the  maximum-dose  relationship  of  Table  2  for  which  b  =  3/2,  a  derivation; 

similar  to  the  one  given  above  for  T  results  in 
&  op 


T  =  0.  5  Z4/5  (px)1/2, 
op  r  '  ' 


T  =1.2  Z4/5  (px)3/5, 

op  r 


Z  <  6 


Z  >  12 


(24a) 


(24b) 


In  general,  the  exponent  of  p  x  in  Equation  (23)  scales  linearly  with  b. 
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The  above  equations  for  T  agree  within  10  percent  to  20  percent  (based  on 
graphical  data)  with  parametric  computer  analysis. 

2.  5  SUMMARY 

The  exponential  function,  Equation  (13),  as  a  general  operative  representa¬ 
tion  of  dose  profiles,  is  presented  in  the  context  of  x-ray  effects  analysis. 

Its  relationship  to  blackbody  spectra  is  described  in  terms  of  the  surface 
dose  defined  by  Equation  (16),  the  characteristic  dose  attenuation  depth  of 
Equation  (17),  and  the  energy  balance  relationship  of  Equation  (20).  It  is 
used  to  derive  the  blackbody  temperature  corresponding  to  the  maximum 
dose  defined  by  Equation  (9),  as  expressed  by  Equations  (23)  and  (24).  When 
dose  profiles  are  needed  for  D/D  ^  0.2  for  specific  blackbody  temperatures, 
parametric  energy-deposition  data  must  be  used  to  obtain  reasonable 
accuracy.  Dose  estimates  for  specific  spectra  can  generally  be  determined 
using  parametric  energy-deposition  data,  as  described  in  Section  2.  2. 

Some  of  the  equations  in  this  and  subsequent  sections  refer  to  a  range  of 
applicable  atomic  numbers.  In  general,  the  atomic  number  refers  to  the 
effective  atomic  number,  Z^.  In  cases  where  a  range  of  is  not  covered 
(e.  g.  ,  6  <  Zg  <  12),  the  situation  is  undefined,  but  often  the  numerical 
difference  between  the  two  equations  is  not  large,  so  that  a  choice  is  not 
difficult  as  to  whether  one  of  the  equations  is  adequate.  If  there  is  insuffi¬ 
cient  reference  to  the  energy- deposition  data  or  a  conservative  rationale  to 
make  one  of  the  two  equations  acceptable,  one  must  resort  to  another 
approach  for  these  transition  situations.  While  for  T  the  two  equations  are 
very  similar  for  the  two  ranges  of  Z^,  so  that  there  is  usually  not  a  problem 
for  6  <  Zg  <12,  Equations  (16)  and  (17)  are  significantly  different  for 
6  <  Z^  <12.  Therefore,  for  this  range  of  effective  atomic  numbers  it  is 
usually  best  to  resort  to  parametric  energy-deposition  data  to  determine 
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Section  3 

FRONT -FACE  MASS  REMOVAL 

Mass  removal  at  the  front  face  of  a  material  may  involve  vaporized, 
melted,  or  virgin  material.  The  mass  removal  plane  is  always  determined 
by  the  phenomenon  of  spall  caused  by  a  thermomechanical  tensile  wave, 
regardless  of  the  degree  of  vaporization  or  melt,  or  the  relative  amplitudes 
of  impulse -induced  stress  and  thermomechanical  stress. 

There  are  two  kinds  of  criteria  for  front-face  mass  removal  which  are 
often  used  to  define  the  spall  plane,  a  dose  criterion  and  a  spall  threshold 
criterion.  The  dose  criterion  corresponds  approximately  to  the  energy 
required  to  reach  the  melt  temperature,  at  which  little  resistance  to  a  ten¬ 
sile  stress  wave  remains.  The  spall  threshold  is  usually  assumed  to  be  the 
dynamic  strength  of  the  unheated  material.  Both  represent  a  simplification 
of  the  real  situation,  and  a  clear  choice  between  the  two  has  not  been  demon¬ 
strated  for  all  situations.  Most  comparisons  have  indicated  that  a  dose  cri¬ 
terion  correlates  satisfactorily  with  experiment  for  many  materials,  although 
there  are  exceptions,  notably  graphite. 

The  modified  BBAY  model  (see  Section  5)  employs  a  dose  criterion 
equal  to  the  energy  through  melt  to  account  for  nonvapor  mass  removal. 

Mass  removal  as  well  as  total  impulse  correlates  reasonably  well  with  this 
criterion  for  situations  involving  significant  vaporization.  Mass  removal 
for  metals  has  been  shown  to  correlate  for  both  x-ray  and  electron  beam 
dose  profiles,  with  a  dose  plane  corresponding  to  between  one  half  the  melt 
temperature  and  the  melt  temperature.  Composite  materials  like  fiber  glass 
have  shown  a  similar  correlation.  This  is  not  inconsistent  with  data  which 
show  that  strength  is  relatively  unaffected  by  temperature  for  very  short 
times  (see  Section  4),  at  least  not  until  the  melt  phase  is  approached.  For 
most  structural  materials,  the  energy  to  melt  and  the  energy  through  melt, 
relative  to  room  temperature,  differ  by  only  a  ratio  of  2/3,  and  often  by  less. 
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However,  graphite,  which  is  unique  in  that  it  has  no  melt  phase  and  its 
strength  increases  with  temperature,  does  not  correlate  well  with  a  dose 
criterion.  On  the  basis  of  existing  experimental  evidence,  a  dose  criterion 
is  recommended,  except  for  graphitic  materials.  Using  a  dose  criterion, 
designated  Dr,  the  maximum  thickness  of  material  removed  is  obtained  from 
Equation  (6),  and  is 


x 


r 


0.  24 


4> 

PDr 


2  3 

where  x^  is  in  cm,  i-n  cal/cm  ,  p  in  g/cm  ,  and  in  cal/g.  The  corre¬ 
sponding  optimum  blackbody  temperature  is  given  by  Equation  (23).  The 
fraction  of  the  incident  fluence  absorbed  in  the  remaining  thickness  of  a 
material  is  approximated  by 


f  = 


(26) 


where  h  is  the  original  material  thickness  in  cm  (see  Sections  2  and  4).  This 
can  be  used  to  estimate  the  average  temperature  rise  of  the  remaining 
material  in  order  to  evaluate  the  strength  reduction  due  to  the  combined 
effects  of  mass  loss  and  temperature  rise  for  situations  where  mass  loss  is 
the  dominant  effect.  For  the  analysis  of  situations  where  temperature  rise 
is  the  dominant  effect,  see  Section  4. 

To  determine  the  mass  removed  for  a  specific  spectrum,  the  use  of 
parametric  energy-deposition  data  is  recommended. 

The  dose  criterion  Dr  is  equal  to  the  energy  required  to  raise  the  target 
enthalpy  to  a  certain  value.  Its  numerical  value  is  therefore  determined  by 
the  ambient  temperature  of  the  target  material. 
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Section  4 

STRENGTH  REDUCTION  BY  HEATING 


One  x-ray  effect  is  the  heating  of  exposed  load-bearing  structures  and 
the  consequent  reduction  in  their  strength  (as  well  as  in  the  strength  of  any 
contiguous  bonds).  There  are  two  major  factors  involved  in  strength  loss 
caused  by  heating:  (1)  the  removal  of  mass  and  the  heat  capacity  and  x-ray 
energy  it  contains  by  front-face  spall,  and  (2)  the  temperature  rise  of  the 
remaining  structure  caused  by  the  residual  x- ray  energy.  The  first  part  of 
the  problem  is  to  determine  the  optimum  conditions  that  will  maximize  the 
strength  reduction  for  a  given  fluence.  The  second  part  is  to  obtain  a  rela¬ 
tionship  for  estimating  the  maximum  strength  reduction. 

4.  1  OPTIMUM  THERMAL  CONDITIONS  FOR  STRENGTH  REDUCTION 
Generally,  strength  reduction  is  a  function  of  the  temperature  profile 
across  the  material  thickness.  The  initial  temperature  profile  created  by 
x-rays  is  exponential-like,  decreasing  monotonically  with  thickness.  After 
a  short  time  (typically  on  the  order  of  seconds  for  metals  and  minutes  for 
composite  materials),  the  profile  becomes  nearly  flat  due  to  thermal  con¬ 
duction.  Therefore,  the  two  extremes  of  temperature  profile  are  the  initial 
profile  corresponding  to  the  x-ray  energy-deposition  profile  and  the  flat 
quasi-equilibrium  profile.  The  question  is  what  gives  the  greatest  strength 
reduction. 


The  three  classes  of  strength-vs.  -temperature  relationships  that  character¬ 
ize  most  structural  materials  are  shown  in  Figure  10,  where  o~r  is  the 
reference  strength  at  or  near  room  temperature,  Tt  is  the  threshold  tem¬ 
perature  where  significant  strength  reduction  sets  in,  and  Tm  is  the  melt 
temperature  approximately  where  the  strength  goes  to  zero. 

The  following  generalizations  can  be  shown  mathematically  for  these 
classes  of  curves.  For  a  convex  curve  (<r"  =  d^cr/dT^  <  q),  the  optimum 
temperature  profile  is  flat,  with  the  temperature  equal  to  T  and  the  tliick- 
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Figure  10.  Typical  Qualitative  Relationships  between  Strength  and  Temperature 


ness  of  the  heated  material  corresponding  to  the  amount  of  energy  absorbed. 
This  is  approximated  by  the  initial  profile.  For  a  concave  curve  (<r"  >  0), 
the  optimum  profile  is  flat  and  extends  over  the  entire  thickness  of  material; 
i.  e.  ,  the  average  or  quasi  -  equilibrium  temperature.  For  a  linear  curve 
( <t "  =  0),  the  strength  reduction  is  independent  of  the  profile  as  long  as  the 
temperature  does  not  exceed  Tm. 

Therefore,  for  convex  strength/temperature  curves,  the  optimum  tem¬ 
perature  profile  is  the  initial  energy-deposition  profile;  otherwise,  it  is  the 
flat  profile  corresponding  to  the  average  temperature  after  thermal 
equilibration. 

A  more  important  consideration,  however,  is  that  conventional  strength/ 
temperature  curves  are  obtained  for  a  heating  duration  of  one  half  hour  while 
the  time  regime  for  x-ray  effects  is  much  shorter,  on  the  order  of  seconds. 
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Experimental  data^  show  that  the  ultimate  tensile  strength  of  2024  aluminum, 
6AL-4V  titanium,  and  301  stainless  steel  all  drop  sharply  after  rapid  heating 
(by  laser),  and  approach  the  conventional  strength/temperature  curve  within 
several  seconds.  As  noted  above,  this  is  also  the  approximate  time  required 
for  the  average  temperature  to  equilibrate  in  typical  structures.  These  data 
imply  that  because  the  short-term  strength  reduction  with  time  is  far  greater 
than  the  effect  of  the  temperature  profile,  strength -reduction  estimates  for 
nuclear  vulnerability  should  be  based  on:  (1)  conventional  strength/ 
temperature  curves,  and  (2)  the  average,  qua  si -equilibrium  or  bulk  tempera¬ 
ture  rather  than  the  initial  temperature  profile. 

Therefore,  it  is  concluded  that  the  optimum  thermal  condition  for  strength 
reduction  corresponds  to  the  maximum  average  bulk  temperature  rise  for  a 
given  incident  fluence.  This  is  also  the  most  vulnerable  situe,tion  for  the 
heating  of  a  bond  at  the  rear  surface  of  an  exposed  material. 

4.  2  OPTIMUM  CONDITIONS  FOR  HEATING  AND  MASS  REMOVAL 
The  above  conditions  consider  only  temperature  effects,  given  a  quantity  of 
absorbed  energy.  However,  optimum  conditions  must  also  cor.  aider  removal 
of  mass  due  to  the  significant  fraction  of  absorbed  energy  that  is  contains  and 
which  is  not  available  for  conducting  into  and  heating  the  remaining  material. 
While  the  loss  of  strength  of  a  structure  due  to  the  amount  of  load -bearing 
material  removed  can  be  significant  for  thin  structures,  for  most  re  all  tic 
situations  the  synergistic  effects  of  strength  reduction  due  to  mass  loss  and 
strength  reduction  due  to  heating  are  not  important.  Thus,  optimum  mass 
removal  and  optimum  heating,  which  occur  under  different  optimum  condi¬ 
tions  (i.  e.  ,  different  x-ray  spectra),  can  be  considered  independently,  as  a 
reasonable  approximation  for  vulnerability  analysis. 

Based  on  the  above  considerations,  the  analysis  of  heating  effects  can  be 
separated  into  two  essentially  independent  situations:  (1)  the  situation  in 
which  the  maximum  mass  is  removed,  where  this  is  the  primary  factor  and 
the  accompanying  temperature  rise  is  of  secondary  importance;  and  (2)  the 


^Interaction  Study,  Vol,  IV  -  Structural  Degradation  by  Short  Time 
Heating,  Air  Force  Weapons  Laboratory,  AFWL-TR-70-]  57,  Vol.  IV 
December  1970. 
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situation  in  which  the  maximum  temperature  rise  is  the  primary  factor  and 
mass  loss,  except  as  it  affects  the  temperature  rise,  is  of  secondary 
importance.  The  latter  is  analyzed  here;  the  other  situation  is  discussed  in 
Section  3. 


4.3  MAXIMUM  TEMPERATURE  RISE 

The  specific  objective  is  to  derive  a  relationship  that  defines  the  maximum 
temperature  rise,  accounting  for  the  energy  and  heat  capacity  lost  by  virtue 
of  the  mass  removed,  with  the  mass  removal  being  determined  by  a  dose 
criterion. 


The  temperature  rise  AT,  in  “C,  is  given  by 


AT 


(27) 


where  4>  (f  -  f  )  is  the  energy  absorbed  in  the  remaining  thickness  of  mate¬ 
rial  (h  -  x^),  c|)  is  the  incident  fluence  in  cal/cm^,  f  is  the  cumulative 
fraction  of  energy  absorbed  to  depth  h,  fr  is  the  cumulative  fraction  of 
energy  absorbed  to  depth  xr,  h  is  the  original  thickness  in  cm,  xr  is  the 
thickness  in  cm  of  material  removed,  p  is  the  density  in  g/cm^,  and  Cp  is 
the  specific  heat  in  cal/g-°C.  Assuming  an  exponential  dose  profile  as  a 
functional  approximation  to  the  actual  profile, 

h 

f  =  1  -  e"  l  (28) 


f  = 
r 


(29) 


and  the  thickness  of  material  removed  is 


x 

r 


(30) 


4-4 


where  i  is  the  characteristic  absorption  depth  in  cm,  Dq  is  the  surface  dose 
in  cal/g,  and  is  the  dose  criterion  in  cal/g  for  mass  removal.  Combining 
Equations  (29)  and  (30), 


pi  D 


Substituting  Equations  (28),  (30),  and  (31)  into  Equation  (27)  gives 


(31) 


phD 
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-  in  y  -  in 
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where  y  =  h /i  has  been  defined  to  form  a  dimensionless  spectral  parameter. 


The  objective  is  to  maximize  AT,  which  is  done  by  setting  the  differential  of 
Equation  (32)  with  respect  to  y  equal  to  zero.  This  gives 


ye-y  (y-1) 


y-fny-^n 


(y-D 
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(33) 


The  optimum  solution  is  y  =  1.  The  other  solution,  ye~Y  =  phDr/cj>,  defines 
the  limiting  conditions  where  x^_  =  h. 

For  y  =  1,  i  =  h,  and  Equation  (30)  becomes  for  the  mass  removal. 


1  <  * 


<  e 


(34) 


where  the  range  of  validity  is  imposed  so  that  0  <  xr  <  h.  Substituting  y  =  1 
into  Equation  (3  2)  gives  for  the  maximum  temperature  rise  AT 
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(35a) 


where  the  range  of  validity  is  taken  to  be  the  same  as  for  Equation  (34).  For 
the  case  of  no  mass  removal,  =  f  =  0  and  Equation  (27)  becomes 

phDr 

=  w;p  =  ^  (»-« -  ~T~)  •  -pnr  s  1  (35b) 

Equations  (34)  and  (35)  represent  the  solution  for  the  maximum  temperature 
rise  and  associated  mass  removal.  They  have  been  compared  to  parametric 
computer  calculations  for  blackbody  spectra  over  the  full  range  of  parameters 
and  found  to  give  the  same  result  to  within  better  than  10  percent  (typically 
5  percent). 

It  is  noted  that  the  optimum  condition  of  £  =  h  corresponds  to  the  case  of 
maximum  dose  at  the  rear  surface  (see  Section  2),  regardless  of  the  amount 
of  mass  lost.  Good  agreement  is  obtained  for  for  blackbody  spectra, 

without  modification  of  the  result  obtained  with  an  exponential  dose  profile 
(as  was  required  for  estimating  maximum  dose),  because  of  the  compensa¬ 
tion  between  the  factors  in  Equation  (35a)  and  the  insensitivity  of  Equa¬ 
tion  (35b)  to  profile  differences. 

Figure  11  shows  a  plot  of  Equation  (35)  expressed  in  dimensionless  terms. 

A  somewhat  simpler  empirical  linear  equation  is  found  from  the  figure  to 
represent  Equation  (35a): 
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Because  the  optimum  temperature  rise  corresponds  to  the  maximum  dose  at 
h,  the  relationships  defining  the  optimum  blackbody  temperature  for  the 
maximum  temperature  rise  are  the  same  as  those  of  Equation  (23)  with  h 
substituted  for  x: 


Top  =  0.  5  Z  (ph)1/3  ,  Z  <  6 

TQp  =  1.  2  Z4/5  (ph)2/5  ,  Z  >  12 


(37a) 

(37b) 
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4.4  SUMMARY 

Conditions  are  defined  which  maximize  strength  reduction  when  both  mass 
removal  and  heating  occur,  but  when  heating  is  the  dominant  effect.  (See 
Section  3  for  situations  when  mass  removal  is  the  dominant  effect.)  Because 
of  the  importance  of  strength  reduction  after  rapid  heating,  the  maximum 
average  temperature  rise  is  the  parameter  that  best  characterizes  the  maxi¬ 
mum  strength  reduction.  It  is  given  by  Equation  (35),  which  is  displayed  in 
Figure  11.  The  optimum  blackbody  temperature  corresponding  to  the  maxi¬ 
mum  temperature  rise  is  given  by  Equation  (37),  and  the  mass  removal  is 
given  by  Equation  (34).  For  specific  blackbody  temperatures,  and  AT  can 
be  obtained  from  energy-deposition  curve  like  those  shown  in  Figures  2 
through  7. 
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Figure  11.  Optimum  Bulk  Temperature  Rise 
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Section  5 

BLOWOFF  IMPULSE 

When  sufficient  radiation  is  rapidly  absorbed  at  the  surface  of  a  material, 
explosive  expansion  of  the  vaporized  mass  creates  an  impulse  load  that  lasts 
a  few  10's  of  nonoseconds  and  can  cause  damage  by  material  response  (e.g.  , 
stress  wave  effects)  or  structural  response.  Two  models  for  predicting 
blowoff  impulse  have  been  studied  and  used  extensively.  Perhaps  still  the 
most  popular  one  is  the  so-called  BBAY  model  developed  by  H.  A.  Bethe  et  al. 
for  the  Air  Force  Special  Weapons  Center  circa  1962.  A  modification  of  the 
BBAY  model  was  introduced  by  S.  L.  Thompson  et  al.  of  Sandia  Laboratories 
in  1968.  The  modified  model  includes  contributions  to  the  impulse  from 
nonvapor  phases.  The  BBAY  and  the  modified  BBAY  (MBBAY)  models  are 
relatively  simple  relationships,  and  both  compare  well  with  results  of  PUFF- 
type  analyses  which  account  for  hydrodynamic  details. 

Exact  calculations  for  a  general  energy-deposition  profile  using  either  model 
are  inconvenient  because  they  require  a  computer  program.  This  section 
presents  generalized  curves  and  explicit  relationships  which  are  convenient 
and  are  adequate  for  practical  analyses.  These  relationships  are  valid  for  a 
large  class  of  energy-deposition  profiles,  including  those  for  blackbody 
spectra. 


5.  1  BLOWOFF  MODELS 
The  BBAY  model  is  given  by 


I 
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1.2 
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(D-E  )  X  dx 
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(38) 
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where  I  is  the  impulse,  xg  the  thickness  of  material  vaporized,  p  is  the 
density,  D  is  the  dose  at  x,  Eg  is  the  effective  energy  of  sublimation  (energy 
through  vaporization),  and  the  factor  1.2  is  an  average  value  of  a  weakly 
varying  parameter  of  the  BBAY  theory. 

For  x-rays,  energy  deposition  is  taken  to  be  instantaneous  so  that  thermal 
conduction  is  neglected.  Although  it  is  not  necessary,  Eg  is  usually  taken  to 
be  a  constant  in  Equation  (38).  However,  since  it  represents  an  enthalpy 
change  caused  by  radiation  absorption,  its  value  is  relative  to  the  initial 
(reference)  temperature  of  the  target  material.  This  can  be  significant  under 
some  conditions.  If  ionization  is  significant,  the  effective  Eg  can  be  greater 
than  the  reference-temperature  Eg;  if  vapor  condensation  is  significant,  it  can 
be  less.  Accounting  for  either  effect  analytically  is  quite  uncertain.  There¬ 
fore,  Eg  is  taken  as  a  constant  and  usually  assumed  to  be  equal  either  to  the 
theoretical  reference-temperature  value  at  1  atmosphere  pressure,  or  if 
available,  an  experimentally  determined  (i.  e.  ,  correlative)  value. 

The  MBBAY  model  is  given  by 


where  E^  is  the  energy  through  melt  and  is  used  as  an  approximation  to  the 
internal  energy  at  the  triple -point  state.  E^  is  influenced  by  the  initial 
target  material  temperature  in  the  same  way  as  Eg,  but  the  influence  is 
greater  because  less  energy  is  involved  in  phase  changes. 

5.  2  GENERALIZED  IMPULSE  RELATIONSHIPS 

The  inconvenience  of  evaluating  Equations  (38)  and  (39)  is  due  to  the  difficulty 
of  handling  analytically  complicated  energy-deposition  profiles  like  those 
associated  with  blackbody  x-ray  spectra.  However,  because  I  is  an  integral 
quantity  which  depends  primarily  on  the  90-percentile  portion  of  the  profile 
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regardless  of  the  fluence  level,  it  is  insensitive  to  the  form  of  the  profile, 
especially  for  thicknesses  for  which  D  is  less  than  10  percent  of  the  surface 
dose.  Consequently,  Equations  (38)  and  (39)  can  be  evaluated  by  assuming 
the  functional  form  of  D  to  be 

x 

D  =  D  e  £  (4  0) 

o 

where  Dq  is  the  surface  dose  and  l  is  the  characteristic  absorption  depth. 
Substituting  Equation  (40)  into  Equation  (38)  and  integrating  gives  for  BBAY 


(41) 


(42) 


and 


(43) 


where  DQe  indicates  the  surface  dose  for  an  exact  exponential  dose  profile. 
Similarly,  Equation  (39)  gives  for  MBBAY 
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In  these  equations,  I  is  in  kilotaps,  4>  is  the  fluence  in  cal /cm  ,  and  E  ,  Em> 
and  D  are  in  cal /g.  E  and  E  are  threshold  doses  designated  generally 
as  E^.  These  equations  are  presented  here  as  general  representations  ox  the 
BBAY  and  MBBAY  models  for  exponential- like  dose  profiles,  and  as  being 
adequate  for  blowoff  impulse  analysis.  Although  the  actual  derivation  of  the 
these  equations  is  not  new,  their  general  validity  and  specific  relationship  to 
blackbody  spectra,  as  described  in  this  section,  is  unique. 


5.  3  APPLICATION  TO  EXPONENTIAL-LIKE  DOSE  PROFILES 
Since  conservation  of  energy  is  important  for  blowoff  impulse  just  as  it  is  for 
maximum  dose,  application  of  the  exponentially  derived  relationships  to 
other  similar  dose  profiles  requires  the  use  of  the  energy  balance  derived  in 
Section  2.  From  an  energy  balance,  Equation  (20)  defines  as 
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P*Do 

4> 


Therefore,  q  is  defined  generally  as 


(47) 


0  = 


P*Do 

4>Et 


(48) 
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where  Dq  and  St  are  parameters  representing  any  particular  dose  profile,  as 

described  in  Section  2.  Consequently,  to  apply  Equations  (41)  and  (44)  to  any 

exponential -like  spectrum  requires  only  a  knowledge  of  the  corresponding 

D  and  f  for  that  spectrum, 
o 

Analytical  expressions  for  Dq  and  St  are  given  by  Equations  (16)  and  (17) 
for  blackbody  spectra.  They  result  in: 

Z  <  6  (49a) 


Z  >  12  (49b) 

5.4  DISCUSSION  OF  EQUATIONS 

Equations  (41)  and  (44)  are  plotted  in  Figure  12  with  q  defined  by  Equa¬ 
tion  (48).  They  are  very  similar  except  for  a  shift  in  the  optimum  value  of 
q,  which  is  6.00  for  BBAY  and  10.2  for  MBBAY.  The  optimum  value,  q  , 
defines  the  optimum  energy  spectrum,  but  as  shown  in  the  figure,  the  optimum 
is  broad  and  relatively  insensitive  to  the  spectrum  (which  may  or  may  not  be 
characterized  by  a  blackbody  temperature). 

From  Equations  (42)  and  (45)  and  Figure  12,  it  is  apparent  that  the  maximum 
impulse  for  a  given  fluence,  corresponding  to  the  optimum  spectrum,  is 
easily  obtained  and  requires  a  knowledge  only  of  (j>  and  E^.,  if  the  optimum 
spectrum  falls  within  the  environment  specifications.  This  is  usually  the 
case  and  represents  most  practical  situations  (except  for  beryllium  and 
carbon).  If  the  optimum  temperature  falls  outside  the  specified  spectral 
range,  only  the  actual  and  (  for  the  limiting  temperature  is  necessary  in 
addition. 

As  far  as  contrasting  the  two  blowoff  models  is  concerned,  one  difference  is 
that  the  optimum  blackbody  temperature  can  be  as  much  as  a  factor  of 


q  =  0.  14 


Z34> 

t3ea 


,  „  z2  d> 

q  “  i  ^  5  j2 

T  '  E, 
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2  greater  for  MBBAY  than  for  BBAY.  The  flat  optimum  indicated  in  Fig¬ 
ure  12  and  its  application  to  a  range  of  spectra  tend  to  make  the  difference 
in  optimum  conditions  of  little  practical  significance.  The  other  difference 
in  the  two  models  is  that  the  maximum  impulse  differs  by 

ZM  /  .  , 

max  _  /  1 . 6 

*B  "  I2’3 

max 

Parametric  computer  calculations  for  blackbody  spectra  indicate  that  Equa¬ 
tion  (50)  is  valid  for  most  materials  (i.  e.  ,  all  the  various  materials  that  were 
checked).  Therefore,  practically  speaking  (i.  e.  ,  within  20  percent),  BBAY 
and  MBBAY  are  interchangeable  except  MBBAY  is  approximately  1.4  times 
higher.  The  exceptions  to  this  will  be  for  materials  for  which  Eo/Em  is 
greater  than  6. 


=  1.4  ±20% 


m 


(50) 
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5.  5  OPTIMUM  CONDITIONS  FOR  MAXIMIZING  IMPULSE 


As  noted  earlier,  general  optimization  of  impulse  is  obtained  from  Figure  12 
in  terms  of  q^  for  a  large  class  of  dose  profiles.  This  section  discusses 
techniques  for  handling  blackbody  spectra  specifically. 

Combining  Equation  (49)  with  q^  =  6  for  BBAY  and  q^  =  10.2  for  MBBAY 
results  in  the  optimum  blackbody  temperature  for  maximum  blowoff.  For 
BBAY, 


TOB  =  °'29Z 


ft) 


1/3 


Z  <  6 


TOB  =  0.6Z  z4/5  A 

\  s 


2/5 


,  Z  >  12 


(51a) 


(51b) 


and  for  MBBAY, 


T 

OM 


0.  24 


Z  <  6 


(52a) 


OM 


=  0.  52  Z 


4/5  (_± 
E 


m 


2/5 


Z  >  12 


(52b) 
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Since  for  a  variety  of  materials  (E  /E  )  =  2  ±  20%,  T_w  —1.5  T^^  for 

s  m'  OM  OB 

Z  £  6  and  T_,  ,  —1.7  T^—  for  Z  -  12, 

OM  OB 


For  a  given  fluence.  Equation  (51)  or  (52)  can  be  used  to  determine  whether 
the  optimum  temperature  is  within  the  range  of  temperatures  to  be  considered. 
If  it  is.  Equation  (42)  or  (45)  is  used.  Otherwise,  Dq  and  f  for  the  limiting 
temperature,  along  with  Equation  (41)  or  (44)  or  Figure  12,  is  used  to  deter¬ 
mine  the  maximum  impulse  under  the  imposed  conditions.  In  fact,  for  any 
specific  nonoptimum  temperature,  the  latter  procedure  is  applicable. 
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5.  6  MAXIMUM  BLOWOFF  IMPULSE  IN  SHIELDED  MATERIALS 
The  maximum  blowoff  impulse  in  a  material  shielded  by  air  or  by  a  non¬ 
contiguous  structure  like  a  shroud  can  be  obtained  by  observing  in  Figure  12 
that  the  maximum  impulse  corresponds  to  the  maximum  q  up  to  q  =  q^  . 

For  q  >  p  ,  this  is  still  a  good  approximation  for  most  realistic  situations. 
The  maximum  r|  is  shown  to  correspond  to  the  maximum  D^,  as  follows. 


The  maximum  dose  is  given  by  Equation  (7),  which  is  set  equal  to  Dq: 


D 

m 


D 

o 


0.  24 

P  x 
r  s  s 


(53) 


where  subscript  s  refers  to  the  shielding  material  and  subscript  a  to  the 
absorber  or  target  material.  In  addition  to  the  integral  relationship  dis¬ 
cussed  in  Section  2,  Equation  (18)  indicates  that  at  depth  x  the  optimum 
effective  linear  attenuation  coefficient  is  equal  to  x  (see  Figures  2  through  7). 
Using  Equations  (18)  and  (53)  gives 


i 


s 


=  X 


s 


0.  244> 
o  D 


(54) 


The  linear  attenuation  coefficients  for  the  shield  and  target  are  related  by 


Ps^  8 
Pa*  a 


* 


So  that  combining  Equations  (54)  and  (55)  gives 


,  _  0.244) 

a  p  D 
r  a  o 


(55) 


(56) 
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Substituting  Equations  (53)  and  (56)  into  Equation  (48)  gives  the  maximum  r) 
for  a  shielded  target: 


=  0.24^2  =  (z3-3/Z3-3)  ,  b=l 

Et  PSXSEt  '  a  8  f 


(57a) 


where  b  is  the  thickness  exponent  of  Table  2,  and  the  appropriate  value  of 
n  =  3.  3  has  been  shown.  A  similar  derivation  gives 


D 

\  =  0-27ir 


0.  lift 

(p  X  3/2Et 
'rs  s  t 


Z2.8/Z2.8 


b  =  3/2 


(57b) 


and  in  general. 


aZ4> 


(p  x  )"b  Zn/Zn 
E^b  vrs  s  \  a  s 


(58) 


where  the  fact  that  H  -  -(dfnD  /dx) 

s  m 


=  x  /b  has  been  used, 
s 


Thus,  Equation  (57)  shows  that  the  maximum  q  is  directly  proportional  to  the 

2 

maximum  Dq.  Therefor e, for  a  given  shielding  mass  of  Pgxg  (g/cm  ), 
Equation  (57)  is  used  to  obtain  q  ,  which  is  used  with  Figure  12  to  derive 
maximum  impulse.  For  the  analytical  model  used  here,  this  optimization  is 
exact  for  q  <  qQ  and  is  a  good  approximation  for  q  >  q  . 

In  Equation  (57),  c|>  is  not  the  fluence  incident  on  the  target,  but  the  fluence 
that  would  be  incident  if  there  were  no  shielding  material.  This  formulation 
makes  it  possible  to  retain  the  concept  of  a  blackbody  source.  The  blackbody 
source  temperature  corresponding  to  Equation  (57)  is  given  by  Equation  (23) 
or  Equation  (24),  depending  on  the  value  of  b. 

5.  7  REVIEW  OF  COMPARISONS  WITH  EXPERIMENT 

Both  blowoff  models  have  been  compared  extensively  with  experimental  data, 
some  published, especially  by  the  Air  Force  Weapons  Laboratory  and 
Sandia  Laboratories,  and  some  unpublished.  The  general  indication  is  that 
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both  models  compare  to  within  a  factor  of  2,  with  the  MBBAY  model 
appearing  to  be  slightly  better.  Since  the  uncertainty  in  the  experimental 
data  is  about  a  factor  of  2,  no  clear  choice  is  indicated. 

On  the  other  hand,  the  BBAY  model  correlates  better  than  the  MBBAY  model 

with  blowoff  impulse  caused  by  high-energy  pulsed  lasers.  Experimental 

laser  blowoff  impulse  data  cover  a  larger  range  of  q  and  characteristic 

absorption  depths  than  do  x-rays  because  thermal  conductivity  becomes  an 

effective  energy  transport  mechanism  as  beam  pulse  widths  increase  from 

picoseconds  to  milliseconds.  Also,  blowoff  data  for  lasers  are  significantly 

more  accurate  and  better  correlated  than  for  x-rays.  The  laser  data  correlate 

well  with  the  BBAY  model  over  a  variety  of  materials  and  for  a  large  range 

of  fluences  and  beam  pulse  widths  (the  effect  of  beam  pulse  width  being 

analogous  to  blackbody  spectrum),  but  the  BBAY  model  gives  results  that 

are  consistently  higher  than  those  obtained  by  experiment  by  a  factor  of 

1.  5  to  2.  Also,  laser  data  indicate  that  the  effective  E  is  slightly  greater 

s 

than  the  room-temperature  value,  while  test  data  indicate  the  opposite  for 
x- rays . 

The  conclusion  is  that  either  model  can  be  justifiably  used,  especially  since 
nothing  better  is  available,  and  the  matter  is  largely  one  of  personal  pref¬ 
erence  since  the  results  they  give  are  quite  similar.  For  the  MBBAY  model, 
however,  E  values  are  less  ambiguous  than  the  effective  E  values  for  the 
BBAY  model,  and  they  are  available  for  more  materials.  Further,  the 
MBBAY  model  gives  more  conservative  estimates  from  a.  system  design 
standpoint.  Also,  the  use  of  Em  by  the  MBBAY  model  for  a  threshold 
energy  is  more  consistent  with  other  data  relating  to  front-face  mass  removal. 

5.8  SUMMARY 

For  almost  any  radiation  source.  Equations  (41),  (42),  (44),  and  (45)  give 
blowoff  impulse  estimates  that  are  more  accurate  than  the  demonstrated 
accuracy  of  the  BBAY  and  MBBAY  models,  as  determined  by  comparison 


5 

R.  W.  Langley,  "An  Analytic  Model  for  Laser  Blowoff  Impulse,  11 
McDonnell  Douglas  Paper  No.  WD  2408,  June  1974. 
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with  parametric  computer  calculations.  Equations  (42)  and  (45)  compare 

with  detailed  computer  calculations  for  blackbody  spectra  generally  to 

within  10  percent,  with  no  more  than  20  percent  difference,  for  fluences 
3  2 

to  greater  than  10  cal/cm  .  Figure  12  has  been  compared  for  non¬ 
optimum  situations  only  for  the  BBAY  model.  For  extremely  high  fluences 
(t)  —  1,000),  the  comparison  is  better  than  30  percent.  For  low  fluences,  the 
agreement  is  similar  except  where  very  little  vaporization  occurs.  How¬ 
ever,  this  regime  is  narrow,  is  inherently  uncertain  in  any  case,  and  is  of 
relatively  little  importance  from  a  practical  standpoint.  Equations  (51)  and 
(52)  compare  with  computer  calculations  to  within  about  20  percent. 

The  relationships  given  here  are  recommended  for  blowoff  impulse  analysis. 

The  equations  are  depicted  in  Figure  12.  The  MBBAY  model  is  slightly 

favored  over  the  BBAY  model  because  it  gives  slightly  larger  impulses  and 

has  a  better  defined  energy  threshold  in  E 

°  m 

The  only  quantities  needed  to  define  the  incident  radiation  source  are  4>,  D^, 

| and  f.  Methods  for  determining  optimum  conditions  for  maximum  impulse 
are  discussed  in  Section  5.  5. 
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Section  6 

STRESS  GENERATION 


Stress  waves  generated  by  radiation  fall  into  two  categories:  (1)  those 
induced  by  blowoff  impulse,  and  (2)  thermomechanical  stress  waves.  The 
two  types  always  occur  to  some  extent  concurrently,  but  for  most  practical 
situations  only  the  maximum  of  one  or  the  other  is  of  interest.  Generally, 
when  impulse -induced  stresses  are  maximized  for  a  given  fluence,  thermo¬ 
mechanical  stresses  are  negligible,  and  conversely  when  thermomechanical 
stresses  are  maximized,  impulse-induced  stresses  are  negligible  —  in  fact, 
they  are  zero.  Which  of  the  two  is  dominant  for  a  given  situation  is  deter¬ 
mined  by  the  radiation  spectrum.  Therefore,  for  the  optimum  conditions, 
the  two  types  of  stresses  and  their  associated  damage  mechanisms  are 
analyzed  independently.  In  some  situations,  combined  stresses  may  be  of 
interest.  Analytical  techniques  for  estimating  stress  generation,  (i.e.,  initial 
stress  conditions  near  the  material  surface  exposed  to  radiation)  are  dis¬ 
cussed  below.  Attenuation  of  stresses  as  they  propagate  through  a  material 
is  discussed  in  Section  7. 


6.  1  IMPULSE -INDUCED  STRESSES 


A  blowoff  impulse,  I,  is  characterized  by  an  average  initial  pressure  or 
stress  amplitude  cr^  with  a  temporal  pulse  width  t,  so  that 


cr 

o 


(59) 


where  <r  is  in  kilobars,  t  is  in  microseconds,  and  I  is  in  kilotaps.  The 

value  of  t  is  determined  either  byithe  energy-deposition  time  or  radiation 

shine  time  t  ,  or  by  r  ,  the  time  required  for  the  vaporized  gas  to  expand, 
s  g 

whichever  is  greater. 
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An  estimate  of  r  is  obtained  from 

g 


T 

g 


(60) 


where  £  is  the  characteristic  radiation  absorption  depth  in  cm  and  C  is  the 
acoustic  velocity  of  the  vapor  in  cm/psec.  The  thickness  of  vaporized  mass 
is  approximated  by  3 £,  which  corresponds  to  the  mass  in  which  most  of  the 
energy  is  absorbed.  The  acoustic  velocity  is  given  by 


where  the  specific  heat  ratio  Y  is  5/3  for  a  monoatomic  gas  and  p  is  the 
density  of  the  vapor.  For  a  pressure  of  cr^,  the  vapor  density  is  the  same 
as  the  solid  density;  i.  e.  ,  p  =  p.  Combining  Equations  (59)  through  (61) 

O 


gives  for  r 


g 


r  =  5400piZ/l  (62) 

g 

3 

where  t  is  in  microseconds,  p  is  in  g/cm  ,  £  is  in  cm,  and  I  is  in  kilotaps. 

g 

The  value  of  £  can  be  obtained  from  an  energy-deposition  profile  (e.g.  s  Fig¬ 
ures  2  to  7)  or  it  can  be  estimated  for  blackbody  spectra  from  Equation  (17). 
For  Z  >  12, 


-4  T  3  ^ 

£  =  4x10  (63) 

3 

where  £  is  in  cm,  T  is  in  kev,  p  is  in  g/cm  ,  and  Z  is  the  material  (effective) 
atomic  number. 
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Substituting  Equation  (63)  into  Equation  (62)  gives,  using  the  same  units. 


T 

g 


9x10 


Z  >  12 


(64) 


In  general,  then,  <y  is  determined  from  Equation  (59)  by  using  the  largest 
value  of  t  corresponding  either  to  the  shine  time  or  to  as  given  by 
Equations  (62)  or  (64). 


For  maximum  blowoff  impulse  for  the  BBAY  blowoff  model,  however,  the 
corresponding  optimum  temperature  for  Z  >  12  is  given  by 


(65) 


where  <j>  is  the  fluence  in  cal/cm  and  Eg  is  the  effective  sublimation  energy 
of  the  material  in  cal/g.  Inserting  Equation  (65)  into  Equation  (64)  gives 
t  for  maximum  blowoff  impulse  for  Z  2  12;  a  similar  derivation  applies  to 
Z  ^  6.  The  results  are: 


go 


-4 

-  2x10 


z12/5 

PI 


Z  >  12 
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g° 


6x10^ 


Z  <  6 


(66a) 


(66b) 


Equation  (66)  shows  that  for  maximum  blowoff  impulse,  t^o  <  0.  02  micro¬ 
second  for  most  practical  situations,  and  thus  t  =  t  for  these  cases.  For 

S  2 

example,  for  I  >  1/2  kilotap,  TgQ  <  0.02  microsecond  for  <70  cal/cm  , 
for  all  materials  (cj>is  lowest  for  tungsten).  For  carbon  and  beryllium, 

Tgo 


0.02  for  <(>  <  10^  cal/cm2. 


The  initial  spatial  pulse  width  6q,  which  is  used  in  conjunction  with  <r  as 
initial  conditions  for  analyzing  stress  wave  attenuation,  is  given  by 
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6  =  2Ct  (67) 

o 

where  C  is  the  acoustic  velocity  of  the  solid.  The  factor  of  2  is  introduced 
in  order  to  express  the  square  wave  represented  by  Equation  (59)  as  a  tri¬ 
angular  wave  form,  which  is  the  equilibrium  shape  to  which  all  initial  pulses 
degenerate  as  they  propagate. 

6.2  THERMOMECHANICAL  STRESSES 

Thermomechanical  stresses  are  caused  by  the  forces  associated  with  the 
thermal  expansion  of  a  material  when  it  is  rapidly  heated  internally  by 
radiation.  This  section  considers  only  situations  involving  no  front-face 
mass  removal.  Because  there  is  no  net  impulse  in  these  situations,  a 
rhermomechanical  stress  wave  always  consists  of  a  compressive  wave  fol¬ 
lowed  by  a  tensile  wave  of  equal  but  opposite  impulse.  The  tensile  wave 
chat  immediately  forms  near  the  surface  of  an  exposed  material  (and  can 
cause  fracture,  spall,  or  debonding)  is  the  main  difference  between  thermo- 
mechanical  stresses  and  compressive  blowoff  stresses,  which  can  form  a 
tensile  stress  only  by  reflecting  at  an  inner  boundary. 

For  instantaneous  energy  deposition,  the  initial  stress  is  compressive  and 
has  the  same  shape  as  the  dose  profile,  and  the  peak  stress  is  determined 
only  by  the  maximum  dose  at  the  surface.  However,  for  a  finite  shine  time, 
the  peak  stress  is  reduced  significantly  below  that  for  an  instantaneous 
exposure  at  the  same  fluence,  and  therefore  shine-time  effects  must  be 
accounted  for. 

6,2.  1  Stresses  for  Instantaneous  Energy  Deposition 

Linear  -elastic  behavior  is  assumed  and  the  dose  profile  is  represented  by 
an  exponential  function.  Conservation  of  energy  is  not  important,  but  an 
accurate  r epresenatation  of  the  initial  portion  of  the  dose  profile  is  of  concern. 
With  the  dose  profile  represented  by 

x 

D  -  DQe  1  (68) 


6-4 


Morland 


(6) 


gives  the  compressive 


stress  to  be 


x 

cr  =  0.042  pFD  e  ^  cosh  ~  ,  x  >  Ct  (69) 

co  o  l 

and  the  tensile  stress  to  be 

_Ct 

cr  =  0.042  pTD  e  1  sinh  ~  ,  x  <  Ct  (70) 

to  r  o  f 

3 

where  0.  042  is  the  conversion  factor  between  cal/cm  and  kilobars,  D  is 

o 

the  surface  dose  in  cal/g,  x  is  the  depth  into  the  material  in  cm,  T  is  the 
dimensionless  Gruneisen  coefficient,  t  is  the  time  after  exposure  in  micro¬ 
seconds,  and  the  stresses  are  in  kilobars. 

At  t  =  0,  a  (0)  =  0  and  cr  (0)  =  0.  042  prD(x).  Equilibrium  stresses  existfor 
to  co 

Ct  Z  31,  for  which  cr  (co)  =  1 /2  c r  (0)  and  o\  (<*>)  is  the  same  shape  and  magni- 

c  o  c  o  to 

tude  as  cr  (co),  but  is  oriented  in  the  opposite  direction, 
co 

The  exact  expression  for  the  positive  impulse  of  the  equilibrium  compressive 
stress  wave,  expressed  in  kilotaps,  is 

/co  .  prD  f  , 

cr  =  0.021  — =  0.021^?  (71) 

co  C  C  C 

o 

where  pfDQ  =  (j)  has  been  used  in  deriving  the  final  expression  because  it 
represents  conservation  of  energy,  and  is  valid  in  this  context.  The  negative 
impulse  of  the  equilibrium  tensile  stress  is  also  given  by  Equation  (71), 
since  the  sum  of  both  impulses  is  zero. 


L.  W.  Morland,  Generation  of  Thermoelastic  Stress  Waves  by  Impulsive 
Electromagnetic  Radiation,  AIAA  Journal,  Vol.  6,  No.  6,  pp  1063-1066, 
June  1968. 
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Although  the  impulse  does  not  change  with  finite  shine  time,  the  peak  equilib¬ 
rium  stresses  are  reduced  considerably  for  typical  values  of  i  because  the 
stress  relaxes  during  the  shine  time. 


6.2.2  Shine -Time  Effect 

For  a  constant  radiation  flux  during  a  shine  timeT,  the  compressive  stress 
for  t  £  t  is  given  by  convolution  of  Equation  (70): 


(T 

c 


O' 

co 


0.  042 


Ct 


x 

~i  .  ,  Ct 

e  sinn  — — 

a 


x  >  Ct 


(72) 


Equation  (72)  shows  that  the  peak  initial  compressive  stress  occurs  at  t  =  t 
and  x  =  Ct,  and  that  the  reduction  factor  for  the  peak  stress  is 


c(t) 

ffco(0) 


(73) 


Equation  (73)  represents  the  effect  of  shine  time  in  reducing  the  peak  initial 
compressive  stress.  This  same  factor  also  applies  to  the  reduction  of  the 
peak  equilibrium  (i.  e.  ,  Ct  =  x»  !  ,  t  »  t)  tensile  and  compressive  stresses 
since  they  are  both  one  half  the  magnitude  of  the  peak  initial  compressive  stress. 

6.2.3  Maximum  Peak  Equilibrium  Stress 

The  peak  equilibrium  stress  that  exists  for  x  £  3 1  is  the  main  quantity  of 
interest,  especially  the  peak:  equilibrium  tensile  stress.  When  a  range  of 
radiation  spectra  is  considered,  the  maximum  peak  stress  for  the  optimum 
spectrum  may  be  desired.  Maximizing  peak  stress  for  a  given  fluence  cor¬ 
responds  to  maximizing  Dq.{,  which  minimizes  the  shine-time  effect  while 
maximizing  the  surface  dose,  but  with  the  constraint  that  front-face  mass 
removal  is  avoided. 

The  peak  equilibrium  tensile  or  compressive  stress  for  instantaneous  energy 
deposition  is,  from  Equations  (69)  and  (70). 

o-  =  0.  021  pTD  (74) 

op  o 
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where  Ct  =  x  »  £ .  Combining  Equations  (73)  and  (74)  gives  the  peak  equilib¬ 
rium  tensile  or  compressive  stress  for  finite  t  : 


or 

P 


0,  01 


Prv 

CT 


( 


(75) 


where  Dg  is  the  dose  at  the  surface  of  the  solid  material,  and  is  equal  either 
to  D  ,  the  surface  dose  of  Equation  (68),  or  to  D  ,  the  mass  removal  criterion, 
whichever  is  smaller.  That  is, 


D 

s 


MIN 


(76) 


Equation  (75)  is  the  general  equation  for  peak  stress  for  any  exponential-like 
dose  profile.  It  compares  closely  with  computer  calculations  for  blackbody 
spectra. 


To  obtain  the  maximum  peak  equilibrium  stress,  Equation  (75)  must  be 
optimized.  As  already  noted,  for  an  exponential  profile,  the  integral  relation¬ 
ship  holds  that  pfDQ  =4>-  However,  parametric  computer  calculations  show 
that  for  blackbody  spectra  there  is  a  maximum  value  for  pfD^ : 


Max  (pDQ£)  =  k<j>  (77) 

where  k  =  0.  05  for  Z  <  6  and  k  =  0.  6  for  Z  »  12.  The  maximum  values  are 
a  function  of  T  and  Z  for  Z  2  12,  with  optimum  T  varying  from  2  kev  for 
aluminum  to  5  kev  for  tungsten.  (While  a  constant  k  is  correct,  it  may  be 
noted  that  it  is  somewhat  inconsistent  with  the  empirical  expressions  for  Dq 
and  2  for  Z  2  1 2.  ) 

Using  a  dose  criterion  D  to  determine  mass  removal,  Equation  (77)  gives 

for  the  optimum  £  consistent  with  no  mass  removal  (i.  e.  ,  D  =  D  =  D  ), 
r  'so  r 


£ 


m 


k<}> 

P°r 


(78) 
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Combining  Equations  (75),  (77),  and  (78)  gives  the  maximum  peak  equilibrium 
tensile  or  compressive  stress,  for  any  fluence,  for  blackbody  spectra: 


<r 

pm 


0.  01 


kr> 

Ct 


(79) 


Equation  (79)  is  plotted  in  Figure  13  in  dimensionless  terms.  The  surface 
dose  always  equals  D  ,  It  is  noted  that  Equation  (79),  except  for  k,  is  insen¬ 
sitive  to  material  variations  since  T,  C,  and  D  are  approximately  the  same 
for  many  materials. 


Combining  Equations  (17)  and  (78)  gives  for  the  optimum  blackbody  tempera¬ 
ture  corresponding  to  Equation  (79): 


T 

op 


1.  5  Z 


1/3 


Z  <  6 


(80a) 

CR79 


Figure  13.  Optimum  Thermomechanical  Stress 
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Z  >  40 


(80b) 


op 


As  indicated,"  Equation  (80b)  agrees  well  with  computer  results  only  for 
Z  >  40.  Figures  3  through  7  can  easily  be  used  to  estimate  T  by  setting 
Dq/<))  =  D^/4)  and  estimating  the  corresponding  T  by  interpolation. 


6.2.4  Initial  Conditions  for  Attenuation  Analysis 

It  is  recommended  that  the  attenuation  of  tensile  stress  amplitude  be  ignored 
for  the  following  reasons;  (1)  there  is  no  adequate  theory  to  handle  it,  (2) 
tensile  stresses  are  necessarily  below  the  relatively  low  spall  threshold  and 
thus  attenuate  weakly  with  distance,  (3)  thermomechanical  tensile  stresses 
are  usually  of  concern  only  for  structures  with  relatively  small  material 
thicknesses  in  which  attenuation  is  limited,  and  (4)  ignoring  the  attenuation 
that  does  occur  results  in  conservative  estimates. 


The  attenuation  of  thermomechanical  compressive  stresses  can  be  analyzed 
in  the  same  manner  as  impulse-induced  compressive  stresses,  as  described 
in  Section  7.  To  do  this  requires  that  cr^  and  6q,  the  initial  peak  (equilibrium) 
compressive  stress  and  the  spatial  pulse  width,  be  properly  specified. 


Section  7  defines  cr 

o 

given  by 


and  6  ,  with  the  former  defined  above  and  the  latter 
o 


6 


o 


2  Cl 
Co 


(81) 


Substituting  Equation  (71)  and  Equation  (72)  with  x  =  Ct  into  Equation  (81) 
gives  the  general  expression  for  as 


6 


o 


2C  t 
2Ct 


(82) 


This  equation  is  used  in  conjunction  with  Equation  (75)  (i.e.  ,  <rQ  =  cr^)  as 
initial  conditions  for  determining  stress  wave  attenuation.  It  compares 


closely  with  PUFF-type  calculations,  which  might  be  expected  since  the 
compressive  stress  is  almost  linear  with  depth  for  0  <  x  <  Ct  and  0  <  t  <  t 
(while  the  tensile  stress  is  zero  for  t  <  t)„ 

For  maximum  peak  stress,  Equation  ( 8 Z )  becomes 


2PD  Ct 
r  r 


,-1 


6  =  2Ct  \l-e 

o 


k<{> 


This  equation  is  used  in  conjunction  with  Equation  (79). 


(83) 


Equations  (82)  and  (83)  differ  from  Equation  (67)  by  the  denominator. 


6.  3  COMBINED  STRESSES 

As  indicated  in  Section  5  (See  Figure  12),  there  is  a  narrow  transition 
between  the  conditions  where  no  blowoff  impulse  occurs  (thermomechanical 
stress  is  dominant)  and  the  conditions  approaching  maximum  blowoff  impulse 
(beyond  which  impulse -induced  stress  is  dominant).  Most  practical  situations 
are  outside  this  transition  regime.  Within  this  regime,  it  is  considered  ade¬ 
quate  for  design  purposes  to  treat  impulse-induced  stresses  and  thermomech¬ 
anical  stresses  separately,  and  to  superimpose  them.  For  the  two  types  of 
concurrent  stresses  to  be  consistent,  should  equal  the  energy  threshold 
of  the  particular  blowoff  impulse  model  used  (i.  e.  ,  =  E^.). 

6.4  SUMMARY 

This  section  presents  analytical  relationships  that  define  the  initial  compres¬ 
sive  stress  and  pulse  width  associated  with  blowoff  impulse;  i.e.,  Equations 
(58)  and  (6  7).  The  temporal  pulse  width  t  is  equal  for  most  situations  to  the 
shine  time,  for  <f><70  cal/cm  .  Otherwise,  it  is  determined  by  the  gasdynamics 
and  is  given  by  Equation  (62),  or  for  optimum  conditions,  by  Equation  (66). 

The  initial  stress  amplitude  is  attenuated  and  the  pulse  width  broadened  as  the 
stress  wave  propagates,  as  discussed  in  Section  7. 
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Equations  (58)  and  (67)  assume  a  square  stress  wave,  which  corresponds  to 
a  constant  radiation  flux.  For  the  situations  where  the  pulse  width  is  deter¬ 
mined  by  the  shine  time,  the  stress  spatial  and  temporal  profile  follows  the 
flux  temporal  profile  (i.e.,  it  has  the  same  time-dependence).  Consequently, 
if  a  nonconstant  flux  profile  is  known,  cr^  can  be  defined  accordingly  by  apply¬ 
ing  Equation  (72),  with  Equations  (58)  and  (67)  still  applicable. 

This  section  also  presents  analytical  relationships  for  the  peak  tensile  and 
compressive  thermomechanical  stresses  formed  near  the  surface  of  an 
exposed  material.  Relationships  are  given  for  a  fixed  dose  profile  and  for 
the  optimum  dose  profile  that  maximizes  the  peak  equilibrium  stress  for 
blackbody  spectra. 

Tensile  stress  wave  attenuation  is  disregarded.  Compressive  stress  wave 
attenuation  should  be  accounted  for  as  described  in  Section  7;  the  relationships 
derived  here  are  used  as  initial  conditions  for  subsequent  compressive  stress 
wave  attenuation. 

The  general  relationship  for  peak  (equilibrium)  tensile  or  compressive  stress 
is  Equation  (74),  and  the  corresponding  initial  pulse  width  is  Equation  (82). 

In  general,  it  is  important  to  impose  the  condition  of  Equation  (76).  Similar 
relationships  for  optimum  spectra  are  Equations  (79)  and  (83),  in  which  the 
condition  of  Dg  =  is  incorporated.  These  equations  include  the  important 
effect  of  stress  relaxation  during  a  finite  shine  time.  One  result  of  the  shine¬ 
time  effect  is  that  peak  thermomechanical  stresses  are  weakly  material- 
dependent,  especially  for  optimum  conditions. 

Impulse -induced  stresses  and  thermomechanical  stresses  occur,  for  the  most 
part,  exclusive  of  each  other;  that  is,  either  one  or  the  other  is  almost  always 
clearly  dominant.  This  is  definitely  the  case  for  all  optimum  situations.  If 
one  wishes  to  treat  them  as  concurrent  but  separate,  as  for  a  specific  black- 
body  temperature,  this  can  be  done  straightforwardly  by  setting  equal  to 
the  energy  threshold  of  the  blowoff  model  used,  and  using  the  appropriate 
fluence  and  H  for  the  material  not  removed  by  blowoff,  as  determined  from 
energy  deposition  curves.  The  impulse -induced  stress  would  then  be  deter¬ 
mined  as  has  been  described. 
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Section  7 

STRESS  WAVE  ATTENUATION  IN  UNDISTENDED  MATERIALS 

The  attenuation  of  the  amplitude  of  a  compressive  stress  wave  as  it  propagates 
into  a  structure  is  one  of  the  main  phenomena  that  must  be  accounted  for  in 
analyzing  stress  wave  effects.  Consequently,  it  has  been  the  subject  of  con¬ 
siderable  experimental  and  theoretical  research  that  has  produced  a  large 
body  of  literature  which  is  heavily  oriented  toward  theory  rather  than  practi¬ 
cality.  A  brief  discussion  is  provided  to  put  the  practical  aspects  of  the 
subject  in  perspective. 

7.  1  BACKGROUND  DISCUSSION 

There  are  two  basic  approaches  to  analyzing  stress  wave  attenuation.  The 
first  one,  the  one  that  has  received  the  most  attention,  involves  the  use  of 
computer  programs  such  as  the  PUFF  code,  which  performs  sophisticated 
calculations  to  account  for  (1)  a  detailed  description  of  the  equation-of- state, 

(2)  an  exact  numerical  solution  of  the  hydrodynamic  equations,  (3)  time- 
dependence  of  imposed  loads,  (4)  temper ature -dependence  of  material  proper¬ 
ties,  (5)  nonidealized  wave  forms,  and  (6)  interference  effects,  as  well  as  other 
factors.  The  second  approach  is  more  analytical  and  is  based  on  the  use  of 
so-called  weak- shock  theory,  which  accounts  for  mechanical  effects  only. 

The  essential  difference  between  the  two  approaches  is  that  the  second  uses 
a  simplified  equation-of- state  involving  only  the  Hugoniot  (i.  e.,  a  nonlinear, 
nonelastic,  hydrostatic  stress- strain  relationship).  It  therefore  cannot 
account  for  the  effects  of  thermodynamics,  porosity,  elastic  strength,  phase 
changes,  strain  rate,  work  hardening,  inhomogeneties,  etc.  Exact  solutions 
of  weak- shock  theory  require  a  computer  code  and  can  therefore  include  most 
other  aspects  of  the  stress  wave  problem.  However,  the  real  advantage  of  the 
second  approach  is  its  application  to  obtain  a  convenient,  explicit,  general 
analytical  solution  to  idealized  stress  wave  problems.  This  approach  makes 
possible  the  practical  application  of  stress  wave  attenuation  theory,  and 
greatly  promotes  physical  understanding.  Neither  of  the  two  approaches  can 
properly  treat  the  attenuation  of  tensile  stress  waves. 
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Although  a  number  of  theoretical  questions  can  be  raised  about  the  validity  of 
PUFF-type  calculations  under  certain  conditions  and  for  certain  kinds  of 
materials,  they  appear  generally  to  agree  with  experiment.  Most  uncertain¬ 
ties  in  results  are  due  to  uncertainties  in  the  equation-of- state.  In  any  case, 
detailed  PUFF-type  calculations  represent  the  only  consistent  basis  for 
assessing  simpler  representations  like  weak-shock  theory. 

A  number  of  comparisons  of  weak- shock  theory  have  been  made  with  PUFF- 

type  calculations  for  various  materials  and  stress  levels,  some  published 

(e.g.,  early  work  by  the  developers  on  the  BBAY  blowoff  model  and  a  1970 

7 

study  by  Bertholf  and  Oliver  )  and  some  unpublished.  General  agreement  is 
good  for  most  materials  and  stress  levels  relevant  to  weapon  systems.  Gen¬ 
erally,  when  the  proper  initial  conditions  are  used,  the  results  for  relatively 
large  thicknesses  and  initial  stresses  (i.  e. ,  thicknesses  greater  than  approxi¬ 
mately  100  times  the  initial  pulse  width  and  initial  stresses  greater  than  several 
hundred  kilobars),  and  for  all  thicknesses  and  low  initial  stresses  (  <  100 
kilobars)  are  approximately  the  same  for  the  two  methods  of  analysis.  For 
smaller  thicknesses  and  large  initial  stresses,  the  deviation  is  relatively 
small,  around  20  to  30  percent  for  most  comparisons  that  have  been 
examined.  An  exception  is  porous  materials,  which  involve  crushing  as  the 
dominant  attenuation  mechanism;  they  must  therefore  be  treated  as  a  separate 
class  of  materials,  and  are  discussed  in  Section  8. 

Thus,  for  most  undistended  or  fully  dense  materials,  the  weak-shock  theory 
is  entirely  adequate  for  practical  application  to  the  analysis  of  stress  wave 
effects.  Problems  that  arise  from  the  approximations  involved  are  more 
often  theoretical  than  practical.  Further,  a  precise  knowledge  of  stress 
waves  is  only  as  useful  as  the  accuracy  with  which  material  damage  criteria 
are  known.  There  are  virtually  no  definitive  stress-wave  damage  criteria 
available  for  most  materials  because  of  the  difficulty  of  measuring  such  data. 
Typical  criteria  that  are  used  for  design  (not  research  analysis)  do  not  even 
explicitly  include  pulse  width  as  a  parameter,  and  peak  stress  criteria  are 
often  not  known  to  better  than  a  factor  of  2  or  3. 

7 

L.D.  Bertholf  and  M.  L.  Oliver,  Approximate  Analytic  Expressions  for  the 
Attenuation  of  a  Triangular  Pressure  Pulse  with  Distance,  Sandia  Labora¬ 
tories,  SC-RR-69-596,  February  1970. 
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A  further  general  indication  of  the  validity  of  weak-shock  theory  for  other 
than  very  high  stresses  where  thermodynamic  effects  occur,  can  also  be 
inferred  from  the  extent  to  which  a  linear  relationship  holds  between  shock 
velocity  and  particle  velocity,  which  is  one  representation  of  the  Hugoniot. 
Departures  from  linearity  can  usually  be  traced  to  porosity,  material 
strength  properties  associated  with  a  high  elastic  yield  strength,  or  phase 
transitions.  At  typical  stresses,  phase  transitions  and  thermodynamic 
effects  are  not  important  in  most  materials.  The  Hugoniot  elastic  limit  is 
seldom  important  because  it  is  operative  only  at  stresses  so  low  that  attenua¬ 
tion  is  too  weak  to  be  very  significant  anyway,  at  least  over  typical  thicknesses 
of  high-strength  materials  (e.g.,  maraged  steel).  And  porosity  is  treated  by 
a  different  technique  than  weak-shock  theory,  as  previously  stated.  Of  the 

g 

74  diverse  materials  included  in  Kohn  ,  a  linear  relationship  is  valid  for 
90  percent  of  them. 

Finally,  when  an  attenuation  theory  like  weak-shock  theory  is  referenced  to 
actual  surface  conditions  at  zero  thickness,  is  asymptotically  correct  or 
approximately  correct  at  large  thicknesses,  and  has  the  correct  fundamental 
behavior,  it  is  difficult  for  it  to  be  in  serious  error  if  used  correctly. 

Therefore,  weak-shock  theory  in  general  and  the  following  general  equation 
derived  from  it  in  particular  can  usually  be  used  with  the  same  confidence, 
from  a  practical  standpoint,  as  PUFF-type  analyses  for  all  but  a  few  mate¬ 
rials.  It  is  recommended  for  general  use. 

7.  2  ASYMPTOTIC  WEAK -SHOCK  THEORY 

Weak-shock  theory,  as  applied  to  x-ray  effects  was  developed  by  W.L.  Bade 
et  al.  for  the  Air  Force  Special  Weapons  Center  circa  1962.  This  represen¬ 
tation  involved  numerical  analysis,  and  was  not  applicable  to  small  thicknesses 
because  the  initial  (surface  loading)  conditions  involved  only  the  impulse 
rather  than  the  initial  stress  and  pulse  width  being  imposed  independently. 

7 

Other  derivations  (e.g.,  Bertholf  and  Oliver  )  treat  the  initial  conditions 
explicitly,  but  also  have  the  disadvantage  of  requiring  numerical  analysis. 

7Ibid. 

g 

Brian  J.  Kohn,  Compilation  of  Hugoniot  Equations  of  State,  Air  Force 
Weapons  Laboratory,  AFWL-TR-69-38,  April  1969. 
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The  following  derivation  differs  from  that  of  Bertholf  and  Oliver  in  that 
the  asymptotic  approximation  of  weak- shock  theory  is  used  here;  that  is,  only 
first-order  terms  are  retained  in  the  relationship  between  stress  and  wave 
velocity.  This  approximation  allows  an  explicit  analytical  expression  to  be 
obtained  for  stress  as  a  function  of  propagation  distance,  and  it  gives  good 
agreement  with  PUFF-type  calculations. 

One  way  of  visualizing  stress  wave  attenuation  is  that  the  front  of  the  wave, 

or  the  loading  wave,  travels  at  a  higher  velocity  than  the  rear  of  the  wave,  or 

unloading  wave.  This  spreads  the  stress  wave  with  time  or  propagation 

distance.  Since  the  area  (impulse)  of  the  stress  wave  must  remain  constant 

to  conserve  momentum  and  if  the  shape  of  the  stress  wave  remains  constant, 

then  the  peak  pressure  will  be  inversely  proportional  to  the  pulse  width.  Any 

initial  stress  wave  quickly  degenerates  into  an  approximately  triangular  pulse 

shape,  with  the  peak  stress  corresponding  to  an  almost  vertical  shock  front 

and  the  stress  decreasing  linearly  to  approximately  zero  at  the  rear  of  the 

7 

wave.  Bertholf  and  Oliver  show  that  deviations  from  an  idealized  triangular 
stress  wave  have  no  significant  effect  on  peak  stress  attenuation. 

For  convenience,  the  S-Hugoniot  mentioned  earlier  will  be  used.  It  has  the 
form 

U  =  C  +  Su  (84) 

where  U  is  the  shock  velocity,  C  is  the  acoustic  velocity  at  zero  stress,  S  is 

the  Hugoniot  slope  for  a  given  material,  and  u  is  the  particle  velocity. 

Relative  to  undisturbed  material,  the  stress  wave  moves  with  velocity  Cg  +  u, 

with  C  given  by 
s 

C2  =  (85) 

B  dp  ' 

where  C  is  the  acoustic  velocity  at  finite  stress  c r  and  density  p  .  The 
s  s 

application  of  the  conservation  of  mass  and  momentum  to  Equation  (84)  gives 

7 

the  stress  as  a  function  of  density  as 
7Op.  cit. 
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Ku 


(86) 


(1  -Sa) 


where  K  is  the  bulk  modulus  (K  =  pC^  -  E/3(l-2v),  with  E  Young's 
and  Poisson's  ratio),  and  ot  is  defined  as 


where  p  is  the  density  at  zero  stress.  Using  Equations  (85)  and  (86) 

Cs  (  J -s2a~) 

C  ‘  ~  <1-S«)2 

Expanding  Equation  (88)  and  retaining  only  terms  linear  in  a  gives 
C 

*  1  +  (2S-1  )a 
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Conservation  of  mass  and  momentum  give  for  u  , 

u  =  cra/p 

From  Equation  (86), 

a 


so  that  combining  Equations  (90)  and  (91)  gives 


and  combining  Equations  (89)  and  (91)  gives 
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Bertholf  and  Oliver,  op.  cit. 


modulus 


(87) 


(88) 


(89) 


(90) 


(91) 


(92) 


I 


C8  =  c  +  C  (2S-1)  -g-  (93) 

Equations  (92)  and  (93)  represent  asymptotic  weak-shock  theory. 

7.  3  DERIVATION  OF  ATTENUATION  EQUATION 

The  front  of  the  equilibrium  (triangular)  stress  wave  moves  with  velocity 
Cg  +  u,  where  Cg  and  u  are  defined  by  the  peak  stress.  The  rear  of  the 
stress  wave  moves  with  velocity  Cg  +  u  =  C,  at  approximately  zero  stress. 
Therefore,  the  stress  wave  broadens  with  time  t  and  distance  x  as 


A6  =  (C  +  u  -  C)  t  =  2StrCt/K  =  2S(rx/K  (94) 

where  6  is  the  pulse  width,  A6  is  the  increase  in  6  in  distance  x,  and  x  =  Ct. 
To  conserve  momentum, 

tr6  =  o-  6  =  2IC  (95) 

o  o 

where  <tq  is  the  initial  peak  stress,  6q  is  the  initial  (triangular)  pulse  width, 
and  I  is  the  impulse  load.  For  x  »6  ,  b~  A6  so  that  Equations  (94)  and  (95) 
give 


(96) 


for  large  x.  Using  Equation  (96),  an  effective  finite  thickness  xq  is  assigned 
to  the  surface  to  introduce  the  initial  stress  condition: 


6  K 
o 

Xo  2Scr 

o 


(97) 


and  x  is  replaced  in  Equation  (96)  by  x  +  xq  to  give 
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1  +  2S 


(98) 


This  is  the  stress  wave  attenuation  relationship  obtained  from  asymptotic 
theory.  The  same  result  is  obtained  with  a  more  rigorous  and  longer 
derivation. 


Equation  (98)  is  for  a  single  material.  It  is  easily  extended  to  handle  a 
laminated  structure  by  following  the  above  procedure  and  assigning  a  new  xq 
at  each  material  boundary.  The  result  after  algebraic  simplification,  is 


q- 

or 

o 


1 


-1/2 


K. 


(99) 


where  the  indicated  sum  is  over  all  material  thicknesses  from  the  surface  to 
the  depth  at  which  the  stress  is  evaluated,  and  T^.  is  the  product  of  the 
individual  acoustic  transmission  factors  for  all  intervening  material 
boundaries.  Tfc  is  defined  by 


Tt  =  n  Ti 

i 


(100) 


and 


-1 


T.  =  2 

l 


1  + 


(pCJj 

<PC>i+l 


(101) 


th 

where  T.  is  for  the  i  boundary  formed  by  laminates  i  and  i+1. 
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The  pulse  width  is  given  by  Equation  (95): 


—  =  —  (102) 

6o  «■ 

8 

Table  3  gives  a  compilation  of  S  values  or  Hugoniot  slopes  taken  from  Kohn  . 

They  do  not  vary  much  with  material,  as  is  the  case  with  the  Gruneisen 

o 

coefficient  r,  because  theoretically  T  =  2S-1  (e.g.,  Glover  ),  which  appears 
in  Equation  (93). 

7.4  COMPARISON  WITH  PUFF-TYPE  RESULTS 

Figure  14  shows,  for  aluminum,  a  comparison  between  asymptotic  weak- 

shock  theory  as  represented  by  Equation  (98),  full  weak-shock  theory,  and 

a  PUFF-type  calculation  made  with  the  WONDY  code  The  last  two  curves 

7 

were  taken  from  Bertholf  and  Oliver  .  The  comparison  is  excellent 
for  all  thicknesses  at  50  kilobars,  with  Equation  (98)  giving,  in  this  case,  the 
same  result  as  WONDY.  The  comparison  is  good  for  all  thicknesses  at 
250  kilobars,  with  Equation  (98)  being  approximately  25  percent  lower  than 
the  WONDY  result  for  2  <  x/  6q  <  100,  and  approaching  the  WONDY  result  for 
other  distance  ratios.  Results  for  1,  000  kilobars  are  shown  only  as  a  matter 
of  interest,  since  stresses  this  high  are  of  practical  significance  only  when 
large  thicknesses  are  involved  (and  for  which  6^  is  typically  of  the  order  of 
10  ^  cm).  Even  at  this  high  initial  stress  level,  Equation  (98)  is  within 
25  percent  of  WONDY  for  x/6q  s  100,  but  falls  to  about  50  percent  of  WONDY 
for  smaller  thicknesses.  Even  the  latter  comparison  is  comparable  to 
uncertainties  in  failure  criteria. 


7 Op.  cit. 

^Op.  cit. 
g 

J.  E.  Glover,  Dynamic  Measurements  of  the  Gruneisen  Coefficients  of 
Tantalum,  Uranium,  and  Uranium -Molybdenum  Alloys,  Sandia  Laboratories 
SC-RR-71-0314,  August  1971. 

*°W.  Herrmann  and  P.  Holzhauser,  WONDY,  A  Computer  Program  for 
Calculating  Problems  of  Motion  in  One  Dimension,  Sandia  Laboratories, 
SC-RR-66-601 ,  February  1967. 
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Table  3 

HUGONIOT  SLOPES 


Material 

S 

Material 

Aluminum 

1.  33 

Gold 

1.  50 

Aluminum  (1060) 

1. 44 

Graphite 

N/A* 

Aluminum 

N/A* 

commercial 

(6061-T6) 

Pyrolytic 

1.  73 

24ST  Aluminum 

1.  34 

,  graphite 

921  T  Aluminum 

1.42 

Hafnium 

1.  12 

Aluminum  (2024) 

1.  37 

Armco  Iron 

3.  14 

Antimony 

1.  54 

Lead 

1.45 

Avcoat  II 

1.  68 

Lucite 

1.  83 

Beryllium 

1.  14 

Magnesium 

1. 24 

Bismuth 

1.  93 

Manganin 

1. 72 

Boron  nitride 

2.  73 

Molybdenum 

1.  24 

Bras  s 

1. 43 

Mylar 

1.63 

Cadmium 

1.  66 

Nickel 

1.45 

Niobium 

1.21 

Boron  carbide 

1.81 

Silicon  carbide 

1.66 

Nylon 

1.  55 

Palladium 

2.  00 

Tungsten 

1.16 

carbide 

Paraffin 

1.53 

Carbon  phenolic 

N/A* 

AVCO  Phenolic 

1.17 

3-D  Carbon 

1.43 

fiber  glass 

phenolic 

GE  Phenolic 

1.  05 

Chromium 

1.  38 

fiber  glass 

Cobalt 

1.  30 

Chopped  nylon 

1.  61 

phenolic 

Copper 

1.49 

T  ape  -  wound 

1.  32 

Durite 

1.41 

nylon  phenolic 

Epoxy 

1.  52 

Quartz  phenolic 

1.  03 

C  -7  Epoxy 

1.  44 

Zinc 

1.  55 

3-D  Quartz 

0.96 

f'henoli  r 

Exon 

* .  bo 

Hi-D  Glass 

N/A* 

Zirconium 

1.  02 

Material 


X  -cut 

Crystalline 

quartz 

Phenolic  refrasil 

Platinum 

Plexiglas 

Polyethylene 

Polystyrene 

Polyurethane 

RAD  58B 

OTWR 

Series  1  24  resin 
Silver 

Stainless  Steel 
Type  304 

Stainless  Steel 
Type  304L 

Steel,  mild  EN3 

T  ant  alum 

T  eflon 

Thallium 

Thorium 

Tin 

Titanium 

Tungsten 

TWSP 

Uranium -3  wt 
pet  Molybdenum 

V  anadium 


*The  linear  relationship  of  Equation  (84)  does  not  hold. 


S 


N/A* 


1. 01 
1.  54 
1. 45 
1.47 
1.  20 
N/A* 
0.  68 
N/A* 
1.89 
1.  54 
1.51 

1. 49 

1.  69 

1.  32 
1.  97 
1.  50 
1.  25 
1.  52 

1.  15 

1.  30 

1.  08 
1.53 

1.  20 


7-9 


CB79 


7.5  SUMMARY 

Based  on  the  reasons  and  comparison  given  above,  and  on  other  similar 
comparisons  with  theory  and  a  few  experiments,  Equations  (98)  through  (102) 
are  recommended  for  general  use.  It  is  expected  that  PUFF-type  analyses 
would  be  required  only  for  occasional  verification  and  for  unusual  or  very 
specific  problems. 

The  equations  are  applicable  to  compressive  stresses  due  either  to  impulse 
loads  or  to  thermomechanical  stresses.  The  definition  of  the  initial  condi¬ 
tions  are  discussed  in  Section  6.  In  all  cases,  the  initial  stress  wave  is 
assumed  to  be  triangular.  The  actual  initial  stress  wave  will  relax  to  the 
equilibrium  triangular  shape  within  approximately  the  time  it  takes  for  the 
differential  velocity  between  the  loading  wave  and  the  unloading  wave  to 
propagate  a  distance  equal  to  the  initial  pulse  width.  For  example,  for  an 
initial  square  wave  of  width  6  use  of  Equation  (94)  shows  that  the  distance 
required  for  the  pulse  to  become  approximately  triangular  is  x^,  given  by 
Equation  (97). 
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Tensile  stress  propagation  cannot  be  analyzed  as  rigorously  as  compressive 
stress  propagation.  Because  tensile  stresses  are  necessarily  limited  to 
amplitudes  less  than  the  fracture  threshold  and  because  the  attenuation  at 
these  levels  is  therefore  relatively  weak,  it  is  recommended  that  attenuation 
of  tensile  stresses  be  ignored  for  practical  analysis.  This  will  tend  to  give 
conservative  estimates. 

The  equations  as  they  stand  apply  only  to  first-pass  or  single -transit  stresses, 
which  are  dominant  for  most  situations.  However,  the  equations  can  be 
straightforwardly  applied  to  estimate  reinforcement  due  to  stress  wave 
reduction.  This  is  relatively  easy  for  first-order  reflections,  but  soon 
becomes  tedious  for  multiple  reflections. 
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Section  8 

STRESS  WAVE  ATTENUATION  IN  DISTENDED  MATERIALS 

One  of  the  methods  for  reducing  blowoff-induced  stress  wave  damage  is  to 
provide  a  porous  or  foam  (i.e.,  distended)  material  that  attenuates  the  stress 
wave  amplitude  to  a  tolerable  level.  The  dominant  attenuation  mechanism 
in  porous  materials  is  the  spreading  of  the  pulse  width  caused  by  their  large 
compressibility,  which  is  far  more  effective  than  the  mechanical  hysteresis 
mechanism  of  undistended  materials. 

Aside  from  the  inconvenience  and  cost  of  a  PUFF-type  analysis  of  stress 
wave  attenuation,  a  paucity  of  experimental  data  on  the  equation-of-state  of 
porous  materials  poses  a  difficulty.  Therefore,  an  analytical  relationship 
that  does  not  require  an  explicit  definition  of  the  equation-of-state  is  highly 
desirable.  Such  a  relationship  is  obtained  using  the  phenomonoiogical  snow¬ 
plow  model  (e.  g. ,  Herrmann^). 

8.  1  THE  SNOWPLOW  MODEL 

The  snowplow  model  assumes  that  the  only  stress  wave  attenuation  mechanism 
is  due  to  the  compaction  of  the  porous  material.  It  assumes  that  the  porous 
material  is  compacted  to  its  solid  density  and  that  no  energy  is  absorbed  by 
the  crushing  action,  but  that  all  energy  from  the  imposed  impulse  load  goes 
into  the  formation  of  a  stress  wave  uniformly  distributed  in  the  solid 
(compacted)  material.  The  yield  or  crush  strength  is  assumed  to  be  negligible 
compared  to  the  applied  pressure.  Thus,  as  the  pressure  pulse  encompasses 
more  compacted  material,  it  broadens  and  the  amplitude  decreases.  Since 
basically  only  the  conservation  of  energy  and  momentum  is  involved,  the 
model  requires  no  explicit  knowledge  of  the  material  equation-of-state. 


W.  Herrmann,  Constitutive  Equations  for  the  Dynamic  Compaction  of 
Ductile  Porous  Materials,  J.  Applied  Phys.  ,  Vol.  40,  No.  6.  ,  pp  2490- 
2499,  May  1969. 
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If  the  material,  as  it  is  compacted  at  a  pressure  p,  is  uniformly  stressed 
to  the  same  pressure,  the  compacted  material  represents  a  square  stress 
wave  with  a  spatial  pulse  width  equal  to  the  thickness  of  compacted  material, 
and  the  impulse  is  equal  to  the  imposed  impulse  load.  With  a  uniform  com¬ 
pacting  pressure  equal  to  the  stress  in  the  compacted  material,  the  compac¬ 
tion  energy  expended  is  equal  to  the  energy  in  the  stressed  solid  material, 
the  latter  being  the  sum  of  the  kinetic  energy  and  the  stored  strain-energy. 
Therefore, 


pxf  = 
r  v 


Pd 


f  )xu 

V 


K 

2 


€2x(l- 


f  ) 

v' 


mu 


(103) 


where  p  is  the  compacting  pressure  acting  over  distance  xf^,  with  x  the  front 
of  the  stress  wave  as  measured  from  the  original  surface,  and  f  the  porosity 
or  void  fraction.  The  mass  of  solid  material  which  the  pressure  pulse 
encompasses  is  m  =  pfl-f^)  x,  where  p  is  the  density  of  the  solid  material, 
u  is  the  particle  velocity,  K  is  the  bulk  modulus  of  the  solid  material,  and 
t  is  the  unit  strain  of  the  solid  material.  The  first  term  on  the  right  side 
of  Equation  (103)  is  the  kinetic  energy  and  the  second  term  is  the  strain 
energy.  The  strain  energy  is  equal  to  the  kinetic  energy  for  a  linear -elastic 
material.  Since  the  stress  <r  is  equal  to  the  pressure  and  the  impulse  I  =  mu, 
Equation  (103)  becomes 


cr 


10 


-3 


pf  (1-f  ) 
r  v  v' 


(104) 


which  gives  the  stress  in  kilobars  at  depth  x  cm,  for  a  porous  material  of 

3 

porosity  f  and  solid  density  p  g/cm  ,  and  an  impulse  of  I  kilotaps.  This  is 
the  stress  as  a  function  of  I  and  x  as  given  by  the  snowplow  model. 

8.  2  STRESS  WAVE  ATTENUATOR  SIZING 

The  it-  i  Jr  ■  J  tliivkness  of  a  rtrcsr  v  :  ve  attenuator  material  is  usually  deter¬ 
mined  by  the  stress  amplitude  transmitted  into  the  underlying  structural 
material.  Rearranging  Equation  (104)  and  letting  h  represent  the  porous 
attenuator  thickness  in  cm. 
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/  (PC)s\l 

h  =  0.  0451 

<r  pf  (1-f  ) 
c r  V  V 

l1  +  <PC>a)_ 

where  the  units  are  as  for  Equation  (104),  cr^  is  the  allowable  stress  criterion 

for  the  structural  material,  pC  is  the  acoustic  impedance,  C  is  the  acoustic 

velocity,  and  the  subscripts  a  and  s  refer  to  the  properties  of  the  solid 

attenuator  material  and  the  underlying  structure,  respectively.  Equation 

(105)  includes  the  stress  wave  transmission  factor  T^_  which  relates  the  stress 

at  the  rear  of  the  attenuator,  or  ,  to  the  stress  transmitted  into  the  adjacent 

structure,  <r  ,  where 
s 


T 
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1  + 


(pC)£ 
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(106) 


The  most  desirable  void  fraction  is  that  corresponding  to  the  least  mass  of 
attenuator  material  required  to  achieve  a  given  stress  reduction.  Although 
no  finite  optimum  f  is  indicated  by  the  snowplow  model,  more  detailed  com¬ 
puter  analyses  indicate  that  it  is  approximately  0.5  for  all  materials. 

8.3  SUMMARY 

The  snowplow  model,  as  represented  by  Equations  (104)  and  (105),  has  been 
compared  with  experiment  and  with  more  elaborate  analyses  (e.  g,  , 
Herrmann^),  and  shown  to  be  qualitatively  correct  and  reasonably  accurate. 
For  transmitted  stresses  greater  than  approximately  10  kilobars,  it  pre¬ 
dicts  for  a  given  attenuator  thickness  stresses  that  are  usually  high  by  less 
than  a  factor  of  2,  with  generally  greater  accuracy  for  higher  stresses.  For 
low  stresses  of  a  few  kilobars,  it  tends  to  predict  stresses  that  are  high  by 
approximately  a  factor  of  2.  This  corresponds  to  predicting,  for  a  given 
stress,  thicknesses  that  are  high  by  40  percent.  It  is  therefore  somewhat 
conservative,  but  adequate  for  preliminary  sizing;  more  detailed  calculations 
are  desirable  for  final  design  only  if  an  adequate  equation-of- state  is 
available. 
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Section  9 

STRUCTURAL  RESPONSE 


Structural  response  refers  here  to  the  dynamic  circumferential  motion  of  a 
structure  responding  to  a  rapidly  applied  transient  load,  usually  resulting 
from  blowoff  impulse.  It  has  potential  importance  only  for  large,  relatively 
fragile  structures.  Thermomechanical  pressure  pulses  due  to  x-ray  heating 
of  liquids  and  their  containers  can  also  cause  structural  response. 

Deformation  associated  with  structural  response  can  in  principle  cause 
debonding,  buckling,  or  rupture.  Solid  propellant  damage  is  not  considered 
a  problem  because  of  the  case-limited  strain  and  deflection  of  typical  high¬ 
stiffness  case  materials.  However,  in  reality  it  is  seldom  a  critical  damage 
mode  because  either  the  loads  of  interest  are  too  small  or  material  response 
damage  modes  are  more  critical.  Material  response  is  usually  critical; 
however,  liquid  tanks  can  be  damaged  by  shock  pres sures  induced  by 
relatively  low  fluences. 

Impulsive  loads  2ire  those  whose  duration  is  short  compared  to  typical  struc¬ 
tural  response  times  (0.  1  to  1  msec).  They  are  induced  by  blowoff  impulse 
or  pressure  pulses  induced  in  liquids.  An  equivalent  impulse  load  can  be 
caused  by  a  thermal  load  due  to  rapid  thermal  expansion  in  the  circumferen¬ 
tial  and  consequently  the  radial  directions.  However,  for  a  given  fluence, 
blowoff  impulse  loads  will  usually  be  more  critical. 

This  section  summarizes  analytical  relationships  suitable  for  general 

12 

analysis.  They  are  not  derived  here,  but  are  taken  from  Langley 


1 2 

R.W.  Langley,  Structural  Re sponse  of  Propulsion  Systems,  McDonnell 
Douglas  Paper  No.  WD  2357.  Presented  at  the  Fourth  Symposium  on 
Nuclear  Survivability  of  Propulsion  and  Ordnance  Systems,  April  1974. 
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9.  1  SHELL  BUCKLING  OR  RUPTURE 

The  general  equation  that  defines  the  critical  impulse  to  cause  shell  failure 
is 


(107) 


where  I  is  the  critical  impulse  in  kilotaps,  h  is  the  load-carrying  shell 

thickness  in  cm,  C  is  its  acoustic  velocity  in  cm/psec,  E  is  Young's  modulus 

in  kilobars,  e  is  the  dimensionless  strain,  «  is  the  strain  criterion  at  which 

c 

failure  occurs,  c  is  the  strain  corresponding  to  any  static  stress  that  may 
exist,  and  <r(e)  is  the  (engineering)  stress  of  the  material  in  kilobars  as  a 
function  of  strain  (i.  e.  ,  the  conventional  stress-strain  curve).  For  buckling, 
corresponds  to  the  str ain-to-yield  stress  and  for  rupture  it  corresponds 
to  the  strain-to-ultimate  stress.  For  a  linear  stres s- strain  curve,  the  term 
in  brackets  becomes  (a-  -  g-  ),  where  the  subscripts  correspond  to  those  for 

the  strains  in  the  equation. 


The  factor  F^  accounts  for  a  multishell  structure,  and  is  given  by 
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where  p  is  density,  subscript  1  refers  to  the  shell  whose  response  is  being 
analyzed,  and  the  other  subscripts  refer  to  other  laminates  making  up  the 
multishell  structure.  Equation  (108)  shows  that  the  response  is  greatest 
(Ic  is  smallest)  when  the  mass  of  the  structure  is  smallest. 

Because  a  liquid  cannot  transmit  a  significant  tensile  stress,  liquid  tankage 
can  be  damaged  by  the  positive  half  of  the  momentum  introduced  by  absorbed 
x-rays,  even  if  no  net  impulse  is  created  by  blowoff.  In  this  case,  the 
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critical  impulse  is  given  by  Equation  (107),  but  the  impulse  is  related  to  the 
fluenco  by  [see  Equation  ('71)J 


I  —  0 .  u2 1 


r, 

f  i  (l  -  f) 

e 


(109) 


where  1  is  _t>  tui&taps,  $  is  in  cal /cm  ,  r  is  the  Gruaeisen  coefficient,  sub¬ 
script  t  refers  tc  the  tank,  material,  subscript  (  refers  to  the  liquid,  and  f  is 
the  fraction  of  the  incident  fluence  absorbed  by  the  tank  wall  (but  does  not 
contribute  to  front-face  mass  removal).  If  blowoff  impulse  occurs,  it  will 
cancel  all  or  part  of  the  impulse  of  Equation  (3.09)  because  the  two  types  of 
impulse  art  oppositely  directed. 


9.2  DE  BONDING 

Debonding  by  structural  response  is  very  unlikely,  even  for  large  fluences; 
material  response  damage  is  usually  far  more  critical.  The  equation  for 
estimating  the  critical  impulse  at  which  debonding  will  occur  in  a  two- 
material  bonded  structure  is 


o-  G(u) 
I  =  - 

C  ...i  u; 


(HO) 


where  I  is 
c 


in  kilobars, 


in  kilotaps,  <rc  is  the  dynamic  bond  strength  or  debond  criterion 
and  Aw  is  the  beat  frequency  defined  as 


Aw 


(111) 


with  the  vibrational  frequency  of  the  two  structural  materials,  designated  by 
the  subscripts,  given  by 


C 


w  — 


(112) 


Where  u  is  in  psec 


and  R  is  the  radius  of  the  cylindrical  structure  in  cm. 


G(w)  is  given  by 
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where  a  is  defined  as 


(113) 
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The  minimum  of  Equation  (110)  is 


I 

c  . 
min 
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(114) 


(115) 


and  occurs  at  a  =  1.  This  indicates  that  the  critical  impulse  is  smallest 
when  the  masses  of  the  two  materials  are  approximately  equal,  although  the 
exact  value  of  the  optimum  value  of  a  depends  on  co^/w^  as  indicated  by 
Equation  (113). 
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DESIGN  EQUATIONS 

This  section  combines  some  of  the  relationships  derived  previously  to  obtain 
equations  in  terras  of  explicit  expressions  for  fluence  thresholds  for  specific 
failure  mechanisms  and  typical  structural  configurations.  Some  equations 
for  dose  and  peak  stress  are  also  summarized.  Other  expressions  can  be 
similarly  obtained.  All  the  threshold  expressions  are  for  optimum  spectra, 
corresponding  to  minimum  fluence  thresholds.  The  equations  are  all  for 
normally  incident  radiation:  Section  10.  16  explains  how  oblique  incidence 
is  handled.  All  Z's  refer  to  effective  atomic  numbers. 

If  the  optimum  blackbody  temperature  for  a  given  effect  is  outside  the  range 
of  specified  temperatures,  then  the  threshold  for  the  appropriate  limiting 
blackbody  temperature  must  be  determined.  For  specific  blackbody  tempera¬ 
tures,  explicit  expressions  for  fluence  thresholds  are  generally  not  possible. 
However,  the  equations,  figures,  and  explanations  given  in  this  report  can 
easily  be  used  to  generate  plots  of  fluence  thresholds  versus  the  values  of 
failure  criteria,  and  the  thresholds  for  specific  criteria  and  blackbody 
spectra  can  thus  be  obtained. 

The  symbols  and  units  are  the  same  as  defined  for  the  identified  constituent 
equations  that  are  combined  to  obtain  a  given  expression.  The  design  equa¬ 
tions,  which  are  listed  in  the  table  of  contents,  are  summarized  in 
Tables  4  and  5.  The  units  used  are  calories,  grams,  centimeters,  micro¬ 
seconds,  kilobars,  kilotaps,  and  kev;  all  fluence  thresholds  are  in  cal/cm^. 

10.  1  MAXIMUM  DOSE  IN  A  SUBSURFACE  MATERIAL 

Equation  (9)  gives  directly  the  fluence  threshold  for  the  maximum  dose  in  a 
subsurface  material: 
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(116) 


4>.  =— ^(px)b  (Zn/Zn) 

Tt  a  r  s  a 

where  Dc  is  the  specified  dose  criterion,  and  the  parameters  a,  b,  and  n 
are  given  in  Table  2.  The  quantity  px  is  the  total  material  areal  density  in 
g/cm^  up  to  the  front  face  of  the  subsurface  material  being  considered,  and 
Zs  and  Za  are  defined  by  Equation  (10). 

The  optimum  blackbody  temperature  is  given  by  Equation  (23)  or  (24), 
depending  on  the  value  of  b. 

10.  2  MAXIMUM  FRONT -FACE  MASS  REMOVAL 

The  fluence  threshold  for  maximum  mass  removal  is  given  directly  by 
Equation  (25): 

4>t  =  4.  2  pDrxc  (117) 


where  xc  is  the  specified  threshold  criterion.  The  optimum  blackbody 
temperature  is  given  by  Equation  (23). 

10.3  STRUCTURAL  HEATING  FOR  OPTIMUM  CONDITIONS 

The  fluence  threshold  for  a  specified  bulk  temperature  rise  can  be  obtained 
from  Equation  (35)  for  optimum  conditions.  However,  because  this  equation 
is  transcendental  in  <|>,  the  fluence  threshold  can  be  most  conveniently 
obtained  from  Figure  11  by  a  trial-and -error  procedure.  The  optimum  black¬ 
body  temperature  is  given  by  Equation  (37). 

10.4  PEAK  COMPRESSIVE  STRESS  IN  A  SUBSURFACE  MATERIAL  DUE  TO 
MAXIMUM  BLOWOFF  IMPULSE 

The  general  equation  for  maximum  blowoff  impulse  is,  from  Figure  12, 


I 

max 


A4> 


(118) 


10-2 


Combining  Equations  (59),  (67),  (99),  and  (118)  gives  the  peak  compressive 
stress  in  a  subsurface  material  due  to  maximum  blowoff  impulse: 
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where  T^_  is  the  overall  transmission  coefficient  defined  by  Equations  (100) 
and  (101),  and  includes  all  boundaries  up  to  the  material  in  which  the  stress 
is  evaluated.  The  symbol  Zh  indicates  the  following  summation  of  param- 
ters  in  Equation  (99): 


2. 
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(120) 


The  summation  is  up  to  the  point  at  which  the  stress  is  evaluated. 


The  optimum  blackbody  temperature  is  given  by  Equation  (51)  or  (52), 
depending  on  the  blowoff  model  used. 


10.  5  BACK-FACE  SPALL  OF  A  SURFACE  MATERIAL  DUE  TO  MAXIMUM 
BLOWOFF  IMPULSE 

Rearranging  Equation  (98)  gives 
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cr 
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(121) 


where  cr  is  the  compressive  stress  at  the  back  face  of  the  surface  material. 
Setting  the  reflected  tensile  stress  cr^  to  the  spall  criterion  <r  , 


(r  =  crT  =  a 
r  r  s 


(122) 
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where  Tf  is  the  reflection  coefficient  between  the  surface  material  and  any 
contiguous  subsurface  material,  and  is  given  by 


(PC)2  -  (PC)1 
Tr  =  (PC)1  +  {PC)2 


(123) 


Subscript  1  refers  to  the  surface  material,  subscript  2  to  the  subsurface 
material,  and  pC  is  the  acoustic  impedance.  If  there  is  no  subsurface 
material,  Tr  =  -1.  If  Tr  <0,  this  indicates  a  reflected  stress  wave  that 
changes  from  compressive  to  tensile,  or  vice  versa. 


Solving  Equation  (121)  for  cr0,  and  substituting  for  a  from  Equation  (122)  and 
for  6q  from  Equation  (67),  <rQ  is  given  by 
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Combining  Equations  (59),  (118),  and  (124)  gives|<)>t  for  causing  back-face 
spall  of  a  surface  material  due  to  maximum  blowoff  impulse; 
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(125) 


The  optimum  blackbody  temperature  is  given  by  Equation  (51)  or  (52), 
depending  on  the  blowoff  model  used. 


10.6  BACK-FACE  SPALL  OF  A  SUBSURFACE  MATERIAL  DUE  TO 
MAXIMUM  BLOWOFF  IMPULSE 

Equation  (99)  generalizes  Equation  (98)  by  letting 
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and 


(127) 


where  Zh  is  defined  by  Equation  (120)  and  T  is  the  overall  translation 
coefficient  defined  by  Equations  (100)  and  (101).  Similarly,  Equation  (125) 
is  generalized  to  give  4>t  for  causing  back-face  spall  of  a  subsurface  material 
due  to  maximum  blowoff  impulse; 
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where  Tj.  includes  all  boundaries  up  to  the  material  in  which  the  stress  is 
being  evaluated;  if  it  is  the  last  material,  Tr  =  -1.  Tr  is  defined  by 
Equation  (123).  Zq  is  summed  up  to  the  stress  plane  of  interest. 


The  optimum  blackbody  temperature  is  given  by  Equation  (51)  or  (52), 
depending  on  the  blowoff  model  used. 


10.  7  BACK-FACE  SPALL  OF  A  SURFACE  MATERIAL  DUE  TO 
MAXIMUM  THERMOMECHANICAL  STRESS 

The  maximum  tensile  stress  at  the  back  face  is  equal  to  the  sum  of  the 
reflected,  attenuated  compressive  stress  and  the  unreflected,  unattenuated 
tensile  stress.  The  former  is  given  by  Equation  (98)  multiplied  by  the 
reflection  coefficient  between  the  surface  material  and  any  contiguous  sub¬ 
surface  material,  with  the  <rQ  of  Equation  (98)  equal  to  the  o-pm  of 
Equation  (79)  and  with  bQ  defined  by  Equation  (83).  The  latter  stress  is 
just  equal  to  c r0  or  <rpm  of  Equation  (79).  Using  these  three  equations  to 
combine  the  tensile  and  reflected  compressive  stresses  gives  the  maximum 
peak  tensile  stress  at  the  back  face: 
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where  Tr  is  defined  by  Equation  (123)  with  its  sign  retained,  and  a  is 
defined  by 


/  2pD  Ct\ 

°  =  1  “  8XP  \-  kf-  j  (130) 

The  fluence  threshold  for  spall  is  obtained  by  plotting  versus  <J>  to  deter¬ 
mine  the  4*  =  <j>t  corresponding  to  o-f-  =  as,  where  4>t  i-s  the  fluence  threshold 
and  (rs  is  the  spall  criterion.  If  there  is  no  subsurface  material,  Tr  =  -1. 

For  spall  to  occur  at  the  back  face  rather  than  near  the  front  face,  the 
amplitude  of  the  tensile  stress,  (rQ,  must  be  less  than  a-g.  If  no  significant 
attenuation  of  the  compressive  stress  occurs  and  Tr  =  -1,  the  upper  limit 
of  the  amplitude  of  the  reflected  compressive  stress  at  the  back  face  is  also 
<r0.  Therefore,  Equation  (129)  is  valid  only  for  <rs  <  or^.  <  2<rs  and  is  conse¬ 
quently  of  limited  importance  relative  to  back-face  spall  from  maximum 
blowoff  impulse.  Nevertheless,  it  is  presented  here  for  the  sake  of 
comparison  with  other  spall  thresholds. 

The  optimum  blackbody  temperature  can  be  obtained  from  Equation  (80)  or 
from  Figures  2  through  7,  depending  on  the  effective  atomic  number. 

10.  8  PEAK  COMPRESSIVE  STRESS  IN  A  NONCONTIGUOUS  SUBSURFACE 
MATERIAL  DUE  TO  MAXIMUM  THERMOMECHANICAL  STRESS 

The  situation  analyzed  here  is  for  a  shielded  (subsurface)  material  in  which 
the  thermomechanical  stress  is  generated.  The  shielded  material  is  non¬ 
contiguous  with  the  shielding  material  (e.g.,  a  component  within  a  shroud). 
The  subsurface  materials  can  generally  form  a  laminated  structure  with  the 
first  material  being  the  material  in  which  the  stress  is  generated. 

It  is  assumed  that:  (1)  there  is  no  front-face  mass  removed  from  the 
shielded  material,  and  (2)  the  maximum  peak  stress  corresponds  to  the 
maximum  surface  dose  in  the  shielded  material.  The  first  assumption  will 
virtually  always  be  true  for  realistic  conditions,  and  the  second  assumption 
is  a  good  approximation  for  all  realistic  conditions. 
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The  peak  equilibrium  stress  is  given  by  Equation  (75).  The  surface  dose  Ds 
is,  by  assumption  (1),  equal  to  D0  of  Equation  (76).  The  maximum  dose  Dm 
is  therefore  given  by  Equation  (53)  (for  b  of  Table  2  equal  to  1),  and  the 
corresponding  l  for  the  shielded  absorber  material  is  given  by  Equa- 
tion  (56).  Equation  (53)  can  be  used  to  verify  that  Dq  <  D^.  Combining 
Equation  (75)  with  Equations  (53),  (55),  and  (56)  gives  the  peak  equilibrium 
(tensile  or  compressive)  stress  generated  in  the  shielded  material: 
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where  the  appropriate  value  of  n  =  3.3  has  been  shown  for  b  =  1.  A  similar 
derivation  gives,  for  b  =  3/2, 


(T 


o 


2.7x10 

Ct(PsXs) 


r<j> 

TTz 


p 

■> 

-  1  -  exp 

3PaCT  fz  2*8\ 

"  p  x  \  a  /  s  / 

r  s  s 

u 

b  =  3/2  (131b) 


and,  in  general, 
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where  P  is  defined  as 
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(3  =  1  -  exp 


(133) 


and  the  general  relationship  between  and  xg  has  been  used: 


1  x 

1  =  -  (dfn  D  /dx)"1  =  (134) 

s  m  b 

The  pulse  width  corresponding  to  the  aQ  of  Equation  (132)  is  given  by 
Equations  (82)  and  (133): 
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Substituting  Equations  (132)  and  (135)  into  Equation  (99)  gives,  for  a  specific 
fluence,  the  maximum  peak  compressive  stress  in  a  generally  laminated 
shielded  structure  due  to  the  (approximately)  maximum  thermomechanical 
stress  formed  in  the  first  shielded  material,  which  is  noncontigous  with  the 
shielded  material: 
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where  psxs  is  the  only  shielding  material  parameter  involved,  is  the 
overall  transmission  coefficient  defined  by  Equations  (100)  and  (101),  and 
is  defined  by  Equation  (120).  Tt  and  include  all  boundaries  and  thick¬ 
nesses  up  to  the  stress  plane. 


The  optimum  blackbody  temperature  is  given  by  Equation  (23)  or  (24), 
depending  on  the  value  of  b. 

10.  9  DEBONDING  OF  A  COMFOSITE/BOND/METAL  STRUCTURE  DUE  TO 
MAXIMUM  BLOWOFF  IMPULSE 

The  tensile  stress  for  this  damage  mode  is  due  to  the  compressive  blowoff 
stress  generated  in  the  composite  material,  transmitted  and  attenuated 
through  the  composite /bond /metal  structure,  reflected  at  the  free  metal 
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surface  as  a  tensile  wave,  and  transmitted  without  attenuation  back  through 
the  metal  and  into  the  bond. 

Equation  (128)  is  modified  to  provide  the  debonding  threshold  for  a  compo¬ 
site/bond/metal  configuration  by  setting  Tr  =  -1  and  Tj.  =  Tjg  ^23  ^32» 
where  T  \2  is  the  transmission  coefficient  between  the  composite  material 
and  the  bond,  T  23  is  the  coefficient  between  the  bond  and  the  metal  sub¬ 
structure,  and  T32  (for  the  reflected  tensile  wave)  is  the  coefficient  between 
the  metal  substructure  and  the  bond.  In  general,  Equation  (101)  indicates 
that 
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Thus,  <t>t  is  given  by 
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where  <rc  is  the  debonding  criterion,  C  and  E^.  are  for  the  composite 
material,  and  E^  is,  from  Equation  (120), 
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The  optimum  blackbody  temperature  is  given  by  Equation  (51)  or  (52), 
depending  on  the  blowoff  model  used. 

10.  10  DEBONDINd  OE  A  COMPOSITE /BOND /METAL  STRUCTURE  DUE 

TO  MAXIMUM  THERMOMECHANICAL  STRESS  GENERATED  IN  THE 
METAL 

The  expression  derived  here  is  for  a  situation  similar  to  that  considered 
in  Section  10.8,  since  the  stress  is  generated  in  the  relatively  higher -atomic  - 
number  subsurface  metal,  and  the  shielding  material  is  the  composite  and 
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the  bond.  Because  the  materials  are  all  contiguous,  the  assumption  of 
Section  10.8  regarding  no  front-face  mass  removal  is  not  relevant. 

The  compressive  stress  generated  in  the  metal  is  transmitted  and  attenuated 
through  the  metal,  reflected  at  the  free  surface,  and  transmitted  without 
attenuation  back  through  the  metal  into  the  bond.  Equation  (136)  is  solved 
for  <j>  in  terms  of  cr ;  is  set  equal  to  T32,  where  T32  (for  the  reflected 

tensile  wave)  is  the  transmission  coefficient  between  the  metal  substructure 
and  the  bond.  In  this  case,  =  S3X3/K2  and  pgxg  =  Plxl  +  P2X2»  where 
subscripts  1,  2,  and  3  refer  to  the  composite,  bond,  and  metal,  respec¬ 
tively.  Thus,  is  given  by 
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where  crc  is  the  debonding  criterion,  T32  is  defined  by  Equatio  n  (137),  and  P 
is  defined  by  Equation  (133). 

The  optimum  blackbody  temperature  is  given  by  Equation  (23)  or  (24), 
depending  on  the  value  of  b. 

10.  11  DEBONDING  OF  A  METAL/BOND/SUBSTRUCTURE  DUE  TO  MAXIMUM 
THERMOMECHANICAL  TENSILE  STRESS  GENERATED  IN  THE 
METAL 

The  maximum  peak  tensile  stress  generated  in  the  metal  is  given  by 
Equation  (79).  Since  tensile  stress  waves  are  assumed  not  to  attenuate,  the 
tensile  stress  in  the  bond  is  given  by  Equation  (79)  multiplied  by  the 
transmission  coefficient,  T  between  the  metal  and  the  bond.  Subscripts  1 
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and  2  refer  to  the  metal  and  bond,  respectively.  T  j2  i-s  defined  by 
Equation  (137).  The  fluence  threshold  $  is  thus  given  by 
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This  is  a  transcendental  equation  for  b-j.  which  must  be  solved  by  successive 
approximation.  However,  the  exponential  term  is  insensitive  to  4>t  and  the 
initial  estimate  of  ignores  it.  The  second  approximation  includes  the 
numerical  value  of  the  exponential  term  which  is  evaluated  using  the  initial 
estimate  of  4>t-  This  process  is  repeated  until  4>t  converges,  usually 
requiring  only  two  successive  approximations;  i.e.  ,  an  initial  estimate  and 
a  second  and  final  estimate. 

The  optimum  blackbody  temperature  is  obtained  from  Equation  (80)  or  by 
reference  to  energy-deposition  curves. 

10.  12  DEBONDING  BY  STRUCTURAL  RESPONSE  TO  MAXIMUM  BLOWOFF 
IMPULSE 

Combining  Equations  (115)  and  (118)  gives  directly  the  minimum  for 
causing  debonding  by  structural  response  to  maximum  blowoff  impulse: 
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The  optimum  blackbody  temperature  is  given  by  Equation  (51)  or  (52), 
depending  on  the  blowoff  impulse  model  used. 

10.  13  SHELL  BUCKLING  OR  RUPTURE  BY  STRUCTURAL  RESPONSE  TO 
MAXIMUM  BLOWOFF  IMPULSE 

Combining  Equations  (107)  and  (118)  gives  directly  b^  for  causing  shell 
buckling  or  rupture  due  to  structural  response  to  maximum  blowoff  impulse: 
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where  the  effective  stress-to-failure  is 
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and  is  equal  to  (crc  -  crst)  for  a  linear  stress-strain  relationship,  where  crst 
refers  to  the  static  stress. 


The  optimum  blackbody  temperature  is  given  by  Equation  (51)  or  (52), 
depending  on  the  blowoff  model  used, 

10.  14  SHELL  BUCKLING  OR  RUPTURE  BY  STRUCTURAL  RESPONSE  TO 
MAXIMUM  NONBLOWOFF  IMPULSE 

Combining  Equations  (107),  (109),  and  (144)  gives  directly  <j>^  for  causing 

shell  buckling  or  rupture  due  to  maximum  nonblowoff  impulse  (i.  e,  ,  no 

front-face  mass  removal  occurs)  in  liquid  propellant  tankage: 
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The  optimum  blackbody  temperature  is  the  temperature  that  maximizes  the 
bracketed  term,  with  the  constraint  that  the  surface  dose  is  less  than  the 
mass  removal  dose  criterion.  There  is  no  equation  for  this,  but  usually  a 
good  approximation  (when  r^/C^.  —  r^/C^  is  simply  to  meet  the  surface  dose 
constraint.  For  this,  the  blackbody  temperature  is  given  either  by  Equa¬ 
tion  (80)  or  Figures  2  through  7,  depending  on  the  atomic  number  of  the  tank 
material. 

10.  15  SIZING  OF  STRESS  WAVE  ATTENUATOR  FOR  MAXIMUM  BLOWOFF 
IMPULSE 

Combining  Equations  (105)  and  (118)  gives  directly  the  thickness  of  a  porous 
stress  wave  attenuator  material  for  maximum  blowoff  impulse: 
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where  is  given  by  Equation  (106).  The  optimum  blackbody  temperature  is 
given  by  Equation  (51)  or  (52),  depending  on  the  blowoff  model  used. 


10.  16  SUMMARY 

The  relationships  for  optimum  spectra  judged  to  be  the  most  useful  are  sum¬ 
marized  in  Tables  4  and  5.  It  is  convenient  in  generalizing  the  analysis  of 
x-ray  effects  to  categorize  materials  as  either  surface  or  subsurface.  In 
this  framework,  Table  4  is  for  surface  materials,  and  Table  5  is  for  sub¬ 
surface  (shielded)  materials. 

As  applicable,  equations  are  shown  for  determining  fluence  thresholds  and 
for  calculating  the  magnitude  of  the  effect  in  question.  For  convenience,  one 
blowoff  model,  MBBAY,  has  been  used.  The  symbol  Z  in  the  equations 
refers  in  general  to  the  effective  atomic  number  Z^.  (See  Section  2.  5  for  a 
discussion  of  situations  involving  6  <  Ze  <  12.  ) 


The  equations,  as  shown  in  the  tables  and  throughout  the  report,  are  all  for 
normally  incident  radiation.  For  oblique  incidence,  the  same  relationships 
hold  exactly  (i.  e.  ,  no  analytical  approximation  is  involved),  except  for 
appropriate  scaling  with  the  cosine  of  the  incidence  angle.  Oblique  incidence 
may  be  a  consideration  for  a  surface  whose  size  is  large  compared  to 
material  thicknesses.  The  equations  are  modified  for  oblique  incidence  as 
follows . 


It  can  be  shown  that  for  all  dose  relationships  that  involve  a  fixed  affected 
thickness,  oblique  incidence  is  treated  exactly  by  the  relationships  for 
normal  incidence  if  the  fluence  is  scaled  as  follows: 


4>o 


d>  COS0 
Tn 


(147) 


where  0  is  the  incidence  angle  as  measured  from  the  normal.  That  is,  for 
normal  incidence,  0  =  0.  The  normally  incident  fluence  in  the  equations  is 
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TABLE  4 

GENERALIZED  DESIGN  EQUATIONS  FOR  SURFACE  MATERIALS 


desiun  equations 

OPTIMUM  BLACKUODY 

equations  , 

REMARKS 

THRESHOLD 

.  EF  F  ECT  ••  • 

1 .  a  Threshold  fur  a  Specified  Peak 

Dose  at  Front  Surface 

DJ3 

0.  = — r  ,z<6 

1  3Z3 

DCT2 

0.  =  ,Z>  12 

60Z 

I .  b  Peak  Dose  at  Front  Surface  for  a 
Specified  Fluencc 

3Z3  0 

DQ=— 6 

60  Z  0 

0o=  ?2  .Z>I2 

Optimum  temperature  is  \ 

lowest  of  interest. 

Peak  dose  is  calculated  for 
specific  temperature. 

Peak  surface  dose  in  any 
surface  material  increases 
with  decreasing 
temperature. 

2.  a  Threshold  for  a  Specified  Amount 
of  Fronl-Facc  Mass  Removal 

0l-4.2DfpX( 

2.  b  Front-Face  Mass  Removal  for 

a  Specified  Fluence  ^ 

Y  _  0.24  0 
pDf 

To  =  0.5Z(pXr)l/3,Z<6 

T0=  1.2  Z4l5(p\)2IS,Z>  12 

3.  a  Threshold  for  a  Specified  Peak 
Equilibrium  Temperature 

Rise 

Sec  Remarks. 

|  3.  b  Peak  Equilibrium  Temperature  ^ 

Rise  for  a  Specified  Fluence 

pXD 

0.37  <  — -  <  1  . 

9 

/  PXDr\ 

.  AT=(pXcp)(-  O' 

pXD 

0  >*■ 

To  =  0.5Z(pX)1/3,Z<  6 

T0-  l.2Z4/5(pX)2/5,Z>l2 

Equilibrium  thermal 
temperature  is  the  average 
temperature  through  the 
material  thickness. 
Equation  is  transcendental 
in  0,  and  AT  vs  0  must  be 
plotted  to  find  0=  0|. 

See  Figure  4-2. 

4.  a  Threshold  for  Back-Face  Spall 

Due  to  Blowoff  Impulse 

O.Jl  SXo  2  Em,/2  I 

0.  - S-23 —  1+  I  1  + 

1  CKT,2  l 

+  (  SXos  )  | 

4.  b  Peak  Back-Face  Reflected  Stress 

Due  to  Blowoff  Impulse  for  a 
Specified  Fluence 

,  6^Tr  /  I.6SX0  V1/2 

"P  ^m,/2\  XCr2Em'/2j 

T=0.24Z(^-)i/3,Z<6 

0 

T  =  0.52  Z4/5(t^-)2/5,  Z  <  12 
m 

r 

Accounts  for  a  bonded 
subsurface  material  via  a 
stress  wave  reflection 
coefficient  Tr.  Negative 

Tr  indicates  tensile  stress. 

5.  a  Threshold  for  Shell  Buckling  or 

Rupture  by  Structural  Response 
to  Blowoff  Impulse 

XF  E  *^2 
„  , ,  Armm 

0t  0.63  »e 

5.  b  Deflection  is  not  a  good  general 
indicator  of  damage. 

Same  as  (4)  above. 

Applicable  for  a  single 
shell  or  multi-shell 
structure.  Shells  must  be 
contiguous  but  not 
necessarily  bonded.  For 
a  linear  stress-strain 
relationship,  oc  =  aQ  — 
osl,  where  o$l  is  the  static 
stress. 

6.  a  Threshold  for  Shell  Buckling  dr 
Rupture  by  Structural  Response 
to  Non-Blowoff  Impulse 

X  f  r,  re  1-1 

V1'0) 

6.  b  Same  as  (5)  above. 

See  Section  10.14. 

Applies  to  liquid  propel¬ 
lant  tanks  only.  Maximum 
non-blowoff  impulse 
implies  no  front-face  mass 
removal.  For  a  linear  stress 
strain  relationship,  oe  = 

°c  ‘  °sr  whcre  nSJ  is  the 
static  stress. 
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TABLE  5 


GENERALIZED  DESIGN  EQUATIONS  FOR  SUBSURFACE  MATERIALS 


DESIGN  l.QUA  HON 


.  rilRI  SIlOLD 


I.  a  Threshold  Cm  a  Specified  Dose* 


i  Thicshold  lor  Back-Face  Spall  ol  a  Bonded 
Subsurface  Material  Due  to  Blowoff 
Impulse  in  >he  Surface  Male  rial 

o.ii  O'hjth,  I 


3.  a  Threshold  Inr  a  Specified  Peak  Compressive 
Stress  in  a  Bonded  Subsurface  Matciial  Due  to 
Blowoff  IrrtpuJse  in  the  Surface  Material 

Same  as  2.  a  with  T  =  I . 


4.  a  Threshold  fur  a  Specified  Peak  Thermo- 
mechanical  Compressive  Stress  in  a 
No n-conliguous  Single  Subsurface  Material, 
or  any  ftrt  or  a  Laminated  Subsurface 
Structure* 


96  b  Cr  <7C  X, 
al*T,  A 

•  I -exp 


-b-l 


afA2  X; 


1/2 


2bpCr  j 

~  1 

k/J 

5.  a  Threshold  lor  Dehonding  of  a  Composite/ 
Bond/Mela)  Structure  Due  to  Blowoff 
Impulse  in  the  Surface  Comi-isiic 
Material 


■r 

f  /2CrT  \  2 

1/2 

rU) 

Tt  _T12T23  JM 


6.  a  Threshold  for  Debonding  or  a  Composite/ 
Bond/Meul  Structure  due  to  Thetmo- 
mcclunical  Stress  Generated  in  the 
Metal  Substiuciurc* 


48bS3X3(p|X|  *P2X2)b  1  od 


.r,K3T, 


,/2WaV 

\  S3X3°dA  / 


1/2 


2hi?3C-,r  /zca  \  1 

PlXl  +f*2X2  \Z”s  /  J 


lTOd  [,  f  2P;IW\T 

,  T.,  CXP  l  Mt  )\ 


1 .  h  Peak  Dose  tor  a  Specified  Fluence* 


>(l"A 

D  =  aCK  XS)“,,^T  J 


I.  h  Peak  Back-Face  Rcllectcd  Stress  m  a  Bonded 
Subsurface  Material  Due  to  Blowoff  Impulse 
in  the  Surface  Material  io<  a  Specified 
Fluence. 


I  60  Tj  Tr  /  t  ,6p  I,  \  - 


■1/2 


OPTIMUM  BLACKBODY 
LOCATIONS 


To-0.5Z(S,)b,1.Z<» 
Tu-  ',z>  12 


To-0.24Z(^-)I/3,Z<  6 
to  =  0.52Z4/s(^b2,s,z>i: 


3.  b  Peak  Compressive  Stress  in  a  Bonded  Sub¬ 
surface  Material  Due  n>  Blowoff  Impulse 
hi  the  Surface  Material  for  a  Specified 
Flucncc 

Same  as  2.  b  with  T  ■  I . 


4.  b  Peak  Thermomechamcal  Compressive  Slrcss 
in  a  Non  -Conliguous  Single  Subsurface  Material, 
or  any  Part  of  a  Laminated,  Subsurface 
Siiucture  ior  a  Specified  Flue  nee* 


a  I’T,  0A 


b  > 


5  b  Peak  Tensile  Stress  in  Bond  of  a  Composite/ 
Bond/Metal  Structure  Due  to  BlowoIT  Impulse 
in  the  Surface  Composite  Material  lor  a 
Specified  Fluence 


I  60T,  /  1 .60  2=  \- 


6  b  Peak  Tensile  Stress  in  Bond  of  a  Composite/ 
Bond/Mcial  Siiucture  Due  <•>  Thermo- 
mechanical  Slrcss  Generated  in  the  Metal 
Substructure  lor  j  Spec  Hied  blue  nee* 


al  ,  I 


3  ’32 


A  0 


7.  a  Threshold  for  Debonding  of  a  Metal/Bond/ 

Substructure  Doe  to  Thermomechanical  Tensile 
Stress  Generated  in  the  Surface  Metal  * 


9<>b(y  <p,X,  ♦u2X2) 


*p;X,|b 


7.  b  Peak  Tensile  Stress  in  Bond  of  a  Mctal/Rond/ 
Substructure  Due  to  Thermomechanical  Tensile 
Stress  Generated  in  the  Surface  Metal  for  a 
Specified  Fluence 


l:,ri  Ti; 

9oC.  r 


I  -  exp 


(-^)l 


Same  at  (I)  above. 


Same  at  (2)  above. 


Same  as  ( I )  above  with 
Z,-P|X,  <p2X2 


See  Section  6.2.3. 


For  any  subsurface  material  shielded 
by  any  number  of  other  materials. 


Back-face  need  riot  be  a  free  surface 
and  is  accounted  for  by  the  siren 
reflection  coefficient. Tr  Negative 
Tr  indicates  icnsile  stress. 


Material  properties  are  for  the  sub¬ 
surface  material  in  which  the  stress 
is  generated,  which  is  the  first  non¬ 
contiguous  shielded  materials. 


Compressive  stress  propagates  Into 
metal ,  reflects  as  ten  site  and  ret  ums 
to  bond.  Subscripts  1 , 2^and  3  refer 
to  composite,  bond,  and  metal  sub¬ 
structure.  respectively. 


Subscripts  1,2,  and  3  refer  to 
composite,  bond,  and  metal  sub¬ 
structure,  respectively.  Shielding 
materials  are  the  compoiite  and  hond. 


Subscripts  I  and  2  refer  to  metal 
and  bond,  respectively.  This  equation 
must  be  evaluated  by  successive 
approximation. 


8.  a  Threshold  for  Pebbnding  by  Structural 
Response  to  BiowoIT  Impulse 

*  l25odtm,/2 
0|.- - — - 


1.  b  Peak  Tensile  Stress  in  Bond  oT a  Composite 
Structure  Due  toSiruclural  Response  to 
Blowoff  Impulse  for  a  Specified  Fluence 

n  M  O.RO  Ai.)  F  — 1/“ 


Same  as  (2)  above. 


Not  likely  to  be  s  critical  failure 
mode. 


*See  Table  2-1  for  values  of  a,  b.  ami  n 
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4>n,  and  4>q  is  the  corresponding  obliquely  incident  fluence.  In  other  words, 

for  a  specified  4>  (the  Effect  column  of  Tables  4  and  5),  the  quantity  to  be 

substituted  for  in  the  equations  is  <f>  cos  0.  Conversely,  for  a  calculated 

threshold  <{>  (the  Threshold  column  of  Tables  4  and  5),  the  normally  incident 

threshold  for  oblique  incidence  is  4^/cos  0.  For  T  ,  substitute  px/cos  0  for 

px  and  E  /cos  0  for  E  .  In  considering  oblique  vs.  normal  incidence,  as  a 
s  s 

general  rule  the  effects  for  oblique  incidence  are  less,  thresholds  are 
greater,  and  optimum  temperatures  for  other  than  blowoff  impulse  are 
greater.  Quantities  that  do  not  change  with  incidence  angle  are:  (1)  the 
surface  dose  for  a  surface  material  because  thickness  is  not  involved,  and 
(2)  the  optimum  temperature  for  blowoff  impulse  because  the  energy /mass 
ratio  must  remain  the  same.  The  fact  that  Tq  for  maximum  blowoff  impulse 
does'  not  change  with  0  is  implied  by  Figure  12,  which  indicates  that  for  a 
fixed  r|,  f  is  proportional  to  4>  regardless  of  0.  By  the  same  token,  for  a 
specific  spectrum  (i.e.  ,  a  fixed  i ),  q  varies  with  cos  0  according  to  how  <f> 
varies  with  cos  0. 
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